UNCLASSIFIED 


AD  NUMBER 


AD819301 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  Aug  1967.  Other  requests  shall 
be  referred  to  AFRPL  [RPPR-STINFO] , 

Edwards  AFB,  CA  93523. 


AUTHORITY 


AFRPL  ltr,  27  Oct  1971 


THIS  PAGE  IS  UNCLASSIFIED 


AFRPIi-TR-67-211  -  VOL  II 


This  Document 
Reproduced  From 
Best  Available  Copy 


PROJECT  SOPHY  -  SOLID  PROPELLANT  HAZARDS  PROGRAM 
Volume  II  --  Apoer.dLxes 


R,  B.  Elweil 
0.  Re  Irwin 
R.  W.  Vail,  Jr. 
Aerojet-General  Corporation 


TECHNICAL  REPORT  AFRPL-TR-67-211  -  VOL  H 


August  1967 


"This  document  is  subject  to  special  export  controls  and 
each  transmittal  to  foreign  governments  or  foreign  nationals 
may  be  made  only  with  prior  approval  of  AFRPL  (RPPR-STINFO) . 
f.  Edwards,  California  93523.” 


Sr'  Air  Force  Rocket  Propulsion  Laboratory 

f  Research  and  Technology  Division 

;  v  #  Edwards  Air  Force  Base,  California 

j  Air  Force  Systems  Command,  United  States  Air  Force 


REPRODUCTION  QUALITY  NOTICE 


This  document  is  the  best  quality  available.  The  copy  furnished 
to  DTTC  contained  pages  that  may  have  the  following  quality 
problems: 

•  Pages  smaller  or  larger  than  normal. 

•  Pages  with  background  color  or  light  colored  printing. 

•  Pages  with  small  type  or  poor  printing;  and  or 

•  Pages  with  continuous  tone  material  or  color 
photographs. 

Due  to  various  output  media  available  these  conditions  may  or 
may  not  cause  poor  legibility  in  the  microfiche  or  hardcopy  output 
you  receive. 


1  1  If  this  block  is  checked,  the  copy  furnished  to  DTIC 

contained  pages  with  color  printing,  that  when  reproduced  in 
Black  and  White,  may  change  detail  of  the  original  copy. 


This  document  contains 
blank  pages  that  were 
not  filmed  ' 


NOTICES 


.(i 


"When  U.  S.  Government  drawings,  specifications,  or  other 
data  are  used  for  any  purpose  other  than  a  definitely  related 
Government  procurement  operation,  the  Government  thereby  in¬ 
curs  no  responsibility  nor  any  obligation  whatsoever,  and  the 
fact  that  the  Government  may  have  formulated,  furnished  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other 
4ata,  is  not  to  b*  regarded  by  implication  or  otherwise,  as  in 
any  manner  licensing  the  holder  or  any  oth«*r  person  or  corpora  - 
.  t;ion,’or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  any  way  be  related 
»■  thereto." 


I 


! 


! 

i 


PROJECT  SOPHY  -  SOLID  PROPELLANT  HAZARDS  PROGRAM 
Volume  II  --  Appendixes 

i  R.  B.  Elwell  .,r: 

j  0.  R.  Irwin  1 

I  R.  W.  Vail,  Jr. 

r 


(■  .  1  ,  ■  ■  ■■  ,  .  . 

"This  document  is  subject  to  special  export  controls  and 
each  transmittal  to  foreign  governments  or  .foreign  nationals 
may  be  made  only  with  prior  approval  of  AFRPL  (RPPR-STINFO) , 
Edwards,  California  93523." 


If 


B*eBBr*>  cm. ow 


TABLE  OF  CONTENTS 

Appendix  No.  Page  No. 

I.  One -Dimens ional  Lattice  Model  of  Deto¬ 
nation  in  Heterogeneous  Explosives . t.  1 

II.  Comments  on  the  RUBY  Code.  . .  25 

III.  Test  CD-96  Data . 39 

IV.  Jetting  Phenomenon. . . 45 

V.  Critical  Geometry  Data,  AAB-3189.  ........ .  6l  ■' 

VI.  Microscopic  Methods  for  Determining 

Void  Content. ............................  *  67 

VII.  Shock  Pressure  Attenuation  in  Plexiglas...  75 

VIII.  Shock  Pressure  and  Wave  Area  Data. ....... .  81 

IX.  Flyer  Plate  Da  :a . 85 

X.  Combined  Blast  Data........... . 91 

DISTRIBUTION. . 137 

DD  FORM  1473 . . . . . . .  .  .  .  147 


X.  Combined  Blast  Data........... . .  91 

DISTRIBUTION. . . . . .  137 

DD  FORM  1473 . . . . . .  .  . .  147 


APPENDIX  I 


One -Dimensional  Lattice  Model  of 
Detonation  in  Heterogeneous  Explosives 

R.  F.  Chaiken 
Technical  Consultant 
Ah,  o  jet-Cericr^l  Coro  oration 
Downey,  California 


INTRODUCTION 

•  _  A; 

I  The  problem  of  defining  the  reaction -zone  in  detonating  hetero¬ 

geneous  explosives  is  of  prime  concern  in  critical  diameter  studies* 
This  arises  from  the  fact  that  the  reaction-zone  length  and  the 
related  detonation  reaction  time  determine  the  space-time  region 
where  chemical  support  of  the  detain  ;i  on  front  occurs.  Hence,  the 
detonation  reaction-zone  and  the  detonation  reaction  time  likewise 
establish,  in  part,  the  ra^es  of  energy  loss  by  side  expansion 
which  would  be  required  to  extinguish  the  detonation. 

,  Eyring  (Ref  erence  I -1)  ,  in  his  classical  grain  burning  theory,  sug¬ 

gested  that  for  granular  explosives  the  detonation  reaction  time 
I  \  can  be  expressed  as 
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where 

Tr  -  =  .  detonation  faction  time  t %* 
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R-  \\  ,p*  average  tadids  of  the  explosive  granules 
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raoe  or  ammonium  nxtrcte  a.ia  ammonium  oercntotate  expJLo- 
liif  burning,  expression  'was  modified  by  Andersen  and 

(Pj'ferencips  ,1-2  ahd.Ifl),  who  suggested  that  >  kr  i|wl  'these. 


^a4es"'i!i«:f MVlih’^iar.  pyrolysis -liatey iji,  of  ;the  solid  under'  high 
surface  heat -flirdcond  it  Ions.  '  ,  f"'  '  1 


]  .ff; T  i$T,V . 


conditions. 


Furthermore,  these  authors  pointed  out  that  the  Ey  ing  grain 
burning  expression  implicitly  assumes  that  i^.iiLiou  of  grain 
burning  occurs  instantaneously  after  passage  of  the  detonation 
front;  however,  under  certain  conditions,  the  time  of  ignition 
may  be  significant  in  establishing  the  detonation  reaction  time. 
In  this  connection,  it  was  suggested  that  the  expression  for 
fr  be  modified  to  include  a  grain  burning  ignition  time,  i.e., 


=  T  +  T 
o  b 


(.  1-2) 


where 


tq  =  time  to  ignition  of  grain  burning 
'V’l.  =  grain  burning  time  (Equation  I-l) 


The  requirement  of  an  ignition  time  in  any  complete  expression 
for  tr>  becomes  clear  when  we  consider  the  case  of  a  heterogeneous 
explosive  in  the  limit  of  Rg  0.  In  thij  limit,  the  explosive 
becomes  homogeneous  in  nature,  and  Tr  from  Equation  1-2  vould 
become  identical  tp  the  ignition  time,  xQ4  Indeed,  this  is  be¬ 
lieved  to  be  the  case  for  liquid  explosives  such  as  nitrome thane 
(References  1-4  and  1-5) . 

However,  while  may  be  clearly  defined  in  the  l^.mit  of  negligi- 
ble  x0  or  negligible  xb,  it  is  certainly  not  clearly  defined 
when  Ta  and  are  of  comparable  magnitude.  In  these  cases ,  the 
simple  addition  rul(e  of  Equation  1-2  does  not  obviously  account 
;  for  the  fact  that  the  ignition  and  grain  burning  reactions  might 
proceed  simultaneously,  rather  than  consecutively;  hence,  its 
applicability  to  defining  the  detonation  reaction-zone  is  ques¬ 
tionable.  . 

It  is  the  purpose  of  this  paper  to  describe  a  model  of  the  deto¬ 
nation  reaction -zone  which  treats  ignition  as  an  integral  Dart 
of  the  detonation  reaction.  The  model  utilizes  the  concepts  of 
boththe  Bowden  hot-spot  ignition  theory  and  the  Eyring  grain 
burning  theory,  and  enables  the  derivation  of  an  expression  for 
Tr  in  which  the  physical  significance  of  xQ  and  its  relationship 
to  Tb  is  more  clearly defined.  The  expression  for  detonation 
reaction  time  also  appears  to  have  direct  applicability  to  the 
detonation  of  solid  composite  propellants. 
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GNb -DIMENSIONAL  LATTICE  MODEL 


Development  of  the.  Model 

In  the  discussions  that  follow,  it  is  considered  that  a  hetero¬ 
geneous  explosive  is  a  lattice  network,  composed  of  initiation 
sites  of  finite  volume  connected  by  a  continuous  medium.  When 
a  detonation  ti-ont  overtakes  a  region  in  the  explosive,  the  in¬ 
itiating  sites  react  to  rorra  hoi  spots.  The  formation  of  r  hot 
spot  is  accompanied  (with  zero  delay)  by  a  reaction  growing  from 
the  hot  spot  and  spreading  into  the  continuum.  Complete  reaction 
of  the  lattice  network  defines  the  detonation  react ion -zone. 

For  a  one-dimensional  treatment  of  the  lattice  model,  the  continuum 
reaction  rate  is  assumed  to  be  described  by  a  constant  linear 
velocity,  B.  The  initiating  sites  may  be  of  different  types  ( i) 
characterized  by  a  length  and  a  reaction  rate  constant  k^. 

■*-he  lattice  at  any  time  t  from,  the  instant  of  passage  of  the  deto¬ 
nation  front  will  consist  of  N^Ct)  sites  of  type  i  per  unit  length 
and  N*x  (t)  -  T  N*^(t)  total  hot  spots  per  unit  length.  Assuming 
hot  spot  formation  to  be  a  first-order  kinetic  process,  the  rate 
of  formation  of  hot  spots  is  given  by 


and 


dN*x(t)/dt  =  kjN^t) 

V-'dN*i(t)  r— i 

dNVt)/dt  = 

i  i 

-1 


( I-3a) 


( I-3b) 


where  k.r  is  the  rate  constant  (sec)“  for  formation  of  the  hot 
spot  ana  (ki)~l  can  be  considered  as  an  ignition  time  for  the  i^'h 
type  of  initiating  site. 

At  time  t  =  ®,  the  number  of  hot  spots  formed  is  simply 

dN*T(  t)  =  j|  dN*,.(®)  =  7  X{N£(  *)'d«  (1-4) 

At  time  t  >♦  ,  the  total  fraction  of  lattice  consumed  by  forma¬ 
tion  of  and  reaction  spreading  from  the  dN*<r(* )  hot  spots  is 

df  (<*>)  =  T  j^2B(  t-p*)  +  Li  J  dN*i(*)  (1-5) 

or  .  > 

df<*)  =  2B(t~*0  7  kiNi(*)d*  ♦  7  L£kiNi(^)d&  (1-6) 


x  lie  TCOtiiJ-  U  ttUtlUil  UJ.  .let  l  L  L.iitii.1*  IB  CuuBUimcu  by  icoCtlOu 

spreading  from  all  hot  spots  formed  for  all  values  of  «%  up  to 
j,  =  t  is  then 


/<b  =  t  r— i  /-*  =  t 

df(  *)  =  2B>  h.  /  (■t-.<r,).N-(*)d* 
*  s  o  ■  XJ  *  -  0  1 

<T)  “  t 

+  y>iki  f  NjCflOdft 


Ni(a)d^ 


(1-7) 
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f(t)  =  2B 


■ijf  o  (t-*)dN*.(^> 


N*  •(  t) 


Eh/  1 


dN*..(*> 


(1-8) 


The  rate  of  lattice  reaction  is  given  by  the  time  derivative  of 
Equation  1-8  which  yields 


-  =  2BN*T(t) 


Lt  dN*i(t) 


(1-9) 


The  detonation  reaction  time,  T^  can  be  defined  as  t  -  xr  when 
f(t)  =  1. 


Integrating  Equation  1-9,  i.e.. 


,f(Tr)  =  1 

f(o)  =  0 


df ( t )  =  2B 


rt  = 
‘/t  -  0 


N*T(  t)dt 


t  =  Tr 


►7k  f  dN*  i(  t ) 

*-r*-  >V  t  =  0 


therefore  yields 


1  =  2B 


Vlt)dt  ♦  ?.  LiN*i(xr) 


(  1-10) 


( 1-11) 


«hich  In  terms  of  a  time -averaged  concentration  of  hot  spots, 
N*,j,  becomes 


2bn*t  LiN*i 


a-i2) 


2BN* 


(1-13) 


where  -t 

r 

“N*t  =-~"  J  N*T(t)dt  (1-14) 

o 

It.  is  interesting  to  note  that  when  the  hot  spots  have  no  dimen¬ 
sion,  and  when  the  average  concentration  of  hot  spots  is  con¬ 
sidered  to  be  the  concentration  cf  explosive  granules.  Equation 
1-13  reduces  directly  to  the  grain -burning  expression  which  was 
suggested  for  ammonium  nitrate  and  ammonium  perchlorate  explo¬ 
sives  (References  1-2  and  1-3),  i.e. , 

Tr  .  *  Rg/3  (1-15) 

The  equivalence  of  Equations  1-13  and  1-15  is  readily  seen  by 
recognizing  that  in  this  case  (N*y)-1  is  simply  the  average  dis¬ 
tance  between  the  centers  of  the  explosive  granules. 

The  question  now  arises  as  to  the  effect  of  finite  ignition  times 
qn  xr.  It  is  apparent  that  a  finite  T0  will  modify  the  value  of 
N*-j(t)  in  Equation  1-14.  Thus  we  must  obtain  an  expression  for 
N*x(t)  which,  from  Equation  1-3,  depends  ui>on  the  time  distribu¬ 
tion  of  initiating  sites. 

Now  the  total  rate  of  decrease  of  type  i  initiating  sites  is 
given' by  %.  ■  " 

-dN .  N.df(t) 

-at  ’  kiNi  ♦  -Ti —  «-«) 

The  first  term  in  the  above  rate  expression  is  simply  the  first- 
order  activation  process  to  form  hot  spots,  while  the  second  term 
is  the  rate  of  destruction  of  initiating  sites  by  action  of  the 


cGiiti.riUuui  react  Lon.  that  is,  when  the  continuum  reaction  over¬ 
takes  an  initiating  site,  that  site  number  is  no  longer  available 
to  form  a  hot  spot. 


Substituting  Equation  1-9  into  1-14  we  obtain 


-dN- 

“dtr-=  kiNi  +  Ni 


r  dN*.(tn 

[2mkr  +Lli  “ dt— J 


( 1-17) 


Summing  over  all  i  and  utilizing  Equation  1-3,  Equation  1-17  be¬ 
comes  a  phaffian  differential  equation  in  more  than  two  variables, 
i.c.  , 


dN, 


+  (1+NtL)  dN*T  +  2RN*,rNTdt  -•  0 


(1-18) 


Where  Nx  ”  7  Ni,  the  total  concentration  of  initiating  sites,  and 
L  =  ?.  dN*i,  the  number  average  length  of  a  hot  soot. 

Unfortunately,  Equation  1-18  does  not  lend  itself  to  ready  solu¬ 
tion.  At  this  time,  however,  an  approximation  approach  will  be 
pursued  which  should  have  some  applicability  to  those  cases  where 
only  sites  with  ignition  times  less  than  Tr  (i.e.,  k£Tr  >  1)  need 
be  considered. 


Approximation  Solution  for  Sites  with  Ignition  Times 
Smaller  than  t he  Detonation  Reaction  Time 


In  the  approximation  approach,  it  is  assumed  that  the  major  loss 
of  initiating  sites  is  through  self  activation  to  form  hot  spots. 
Hence,  Hquatron  1-16  can  be  written  as 


-dN. 

_ _ i_ 

dt 


SK  k^N^ 


or 


‘/There  N 
site 


exp  ( -k^t) 


Cl -19) 
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for  N*x>(  t)  yields 


N*t<  t )  = 


V t>  =  Nio 

ia  the  inicial  concentration  of  the  ith  type  initiating 
stitution  of  Equation  1-20  into  Equation  1-4  And  solving 


■4c  •  t 


(I -21) 
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The  average  concentration  of  hot  spots  is  given  by 


N*t  =  (lArr) 


N  •  T  + 
to  ~~ 


r— i  N- 

Tt7 


exp  ( -kx  Tr)  -1 


(1-22) 


Comparing  the  above  expression  with  Equation  1-13  finally  yields 


__1 
2  BN 


To 
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r — 'F  -  r 

“I 

1  -  T.  V^V 

^  [> 

.  -  exp  ( -kj'/V 

(  1-23) 

where  N<j>0  -  7,(N3-)0,  and  =  (N^)0/N>p0  is  the  initial  fraction  of 
the  ith  type  of  initiating  site.  It  should  be  emphasized  that 
the  only  sites  to  be  considered  in  Equation  1-23  are  those  for 
which  k^'T:r  >  1,  since  it  is  only  for  these  sites  that  the  approxi¬ 
mation  of  Eauation  1-19  can  be  valid.  Within  this  constraint,  the 
value  of  N*£  ('tr-)  is  ^  (N£)0  (i-e.  ,  each  site  leads  to  a  hot 
spot),  and  the  term  [l  -  exp  ( -k^Tr) J  is  *  1.  Hence,  a  reason¬ 

able  approximation  to  Equation  1-23  is 


T 


r 


LiNio  I  +  V 
NTo  J  ^  ^To 


(1-24) 


or 


+ 


(1-25) 


Here,  n  is. the  number  of  different  types  of  initiating  sites  and 
S  =  niNx0,  L  =  ?LiNi0/Nx0»  t-0  =  7  Ni0/^iNT9  are  simply  number- 
weighted  average  parameters  of  the  initiating  sites;  viz. ^  the 
distance  hetveen  centers  (S),  the  Length  (L),  and  the  ignition 
time  (t0). 

The  above  expression  for  the  detonation  reaction  time  is  now  in  a 
for  i  that  can  be  readily  applied  to  various  explosive  situations. 


APPLICATION  OF  APPROXIMATION  APPROACH  TO  EXPLOSIVE 
SYSTEMS  CONTAINING  ONE  AND  TWO  TYpES  OF  INITIATING  SITES 


Case  of  a  Single  Type  of  Initiating  Site  ( n  ~I ) 

Let  us  consider  tn  expl  -ivs  with  only  one  type  of  initiating  site 
which  meets  the  criterior  >  1.  In  this  esse  n  *  1  and  the 

number  weighted  average  site  parameters  in  Equation  1-25  become 


r.~u’ 

V  ff  J  r 


'flu*  &&&* 


tlllU  LC><-  >3  > 


2B 


( 1-26) 


The  similarity  between  the  above  expression  and  Equation  1-2  is 
obvious,  suggesting  that  th>  simple  addition  rule  is  reasonably 
applicable  in  this  case.  it  should  be  noted  that  Equation  1-26 
extrapolates  to  the  expected  ar  for  both  negligible  ignition  time 
and  for  negligible  grain  burning  time.  The  latter  exrapolation 
is  expected  to  be  applicable  to  the  case  of  homogeneous  explosives, 
which  can  now  be  considered  as  those  explosi/es  where  VS  -  L)  »  0. 
This  might  be  expected  when  S,  and  hence  L,  approach  molecular 
dimensions. 


In  the  case  of  finite  size  initiating  sites  (  t0  still  negligible) 
Equation  1-26  is  identical  to  the  detonation  reaction  time  ex¬ 
pression  used  at  Aerojet  (Reference  1-7)  to  describe  deconation 
in  2  to  107o  RDX-adulterated  solid  composite  propellant.  In  this 
case,  it  was  considered  that  the  RDX  particles  would  act  as  in¬ 
itiating  sites,  and  that  the  propellant  medium  separating  the  RDX 
would  act  as  the  grain  burning  continuum.  The  reaction  time  ex¬ 
pression  was  derived  in  terms  of  the  mass  fraction,  fRDX»  of 
RDX  particles,  i.e., 

1/3  1 


"'RDX 

2B 


L  V6P  fRDx/ 


-1 


(1-27) 


where  and  Prdx  are  the  propellant  and  RDX  densities,  respec¬ 
tively,  and  LroX  the  diameter  of  the  RDX  particles.  It  can  be 
readily  shown  that  for  spherical  RDX  particles  the  bracketed  term 
in  the  above  expression  is  identical  with  ( S/L  -  1)  . 

It  is  noteworthy  that  the  use  of  Equation  1-27  has  led  to  a  suc¬ 
cessful  correlation  of  the  critical  diameter  of  solid  composite 
propellants  as  a  function  of  f rdx  in  the  range  0.02  <  f  s?  0.1 
(Reference  1-7). 


ASince  N-pc  is  the  number  of  RDX  particles  per  unit  length,  it  is 
also  equal  to  the  cube  root  of  the  number  of  RDX  particles  per 
u  lit  of  propellant  volume. 


'!  ■**-" 


Cqwp  of  Two  Typps  of  Initiating  Sites  (n  =  2) 

Foi  this  case,  Equation  1-25  can  be  expressed  directly  as 


T  = 


r  -  23  (Nlo  +  N2o) 


Fi(Ll  -  l2)  +  L2 


+  Fi  <Tlc  “  t2o^  +  t2o 


where  F£  -  Nio/N-rQ  is  me  nunuaer 
tially  present  (i.e.,  ?,  F^  =  1). 


(1-28) 

is  the  number  fraction  of  sites  of  type  i  ini- 


It  is  interesting  to  examine  the  above  expression  in  light  of  the 
previous  work  at  Aerojet  (Reference  1-7),  in  which  Equation  1-27 
was  modified  to  account  fer  the  effect  of  inherent  initiating  sites 
on  the  detonation  reaction  time  of  solid  composite  propellant 
containing  only  small  amounts  of  RDX  (i.e.,  fRDX  <  0.02).  Here, 
it  was  suggested  that  the  inherent  sites  (e.g.,  flaws,  voids) 
acted  as  an  additional  weight  fraction  c  of  RDX  particles.  There¬ 
fore,  Equation  1-27  becomes 


.  lrbx 

2B 


/ _ — -P.SPX _ 

\  69  (fRDX  +  c)  / 


1/3 


or 


(fRDX  +  c) 


1/3 


-  C, 


(1-29) 


(1-30) 


In  the  same  terminology,  Equation  1-28  (assuming  negligible  igni¬ 
tion  times)  yields 

htUX.  r  PRDX/6  P  )  1/3  fRDX1/3  *  c’  <LX/Lanx)  1  T  ... 
Tr — W  >  -  a/3  .  _. - 72 — 7.  I  a'31) 


t 


J 


where 
tho 


e  c'  =  (fx  pRDx/px)1'  3  ( Lqpy/L^) t  the  subscript  x  referring  to 
inherent  initiating  sites.  Since  Equation  1-31  has  th,e  form  of 

p7T7 


<«X/3rmc  *  V 


(f 


(1-32) 


RDX' 


“MAKS#? 


It  is  evident  that  Equation  I-3L  and  the  Aerojet  expression  given 
by  Equations  I--29  or  jl-3u  are  not  equivalent.  U  the  inherent 
sites  in  Equation  I-3I  are  assumed  to  have  a  diameter  and  density 
equivalent  to  that  of  RDX ,  the  1-D  lattice  expression  for  xr  re¬ 
duces  to 


Tr  =' 


-  ^DX  (  ^  P  RDX7"6  P  ) 


(  f  RDX'' 


-  1  = 


(fRDX5 


(1-33) 


This  expression  while  similar  to  Equation  1-30  is  still  not  iden¬ 
tical  to  it  except  in  the  limit  of  fRjjx  -  0.  A  more  detailed 
quantitative  comparison  between  these  various  reaction  time  ex¬ 
pressions  will  be  given  in  a  later  section  of  this  report  when 
critical  diameter  data  are  discussed. 


It  was  emphasized  earlier  that  the  general  Equations  1-24  and 
1-25  were  applicable  only  to  systems  with  initiating  sites  that 
form  hot-spots  in  a  time  less  than  the  reaction  time  (i.e., 
k^TTj-.  >  1).  Also,  it  was  shown  in  the  previous  sections  that  the 
reaction  time  will  generally  increase  as  the  concentration  of  one 
type  of  site  decreases  (e.g.,  the  case  of  RDX  adulterated  propel¬ 
lant).  Hence,  a  possible  situation  presents  itself  in  which  an 
explosive  with  many  potential  types  of  initiating  sites  has  only 
one  type  of  site  active  at  short  reaction  times,  but  as  conditions 
change  to  increase  the  reaction  time  (e.g.,  decrease  in  loading 
density  or  decrease  in  prime  hot-spot  concentration),  other  sites 
become  active.  This  phenomena  may  have  implications  with  respect 
to  the  contrasting  behavioral  patterns  in  high  explosives  that  is 
discussed  in  Price  (Reference  1-8)  as  well  as  with  respect  to  the 
role  of  inherent  impurities  in  th*s  case  of  RDX -adulterated  pro¬ 
pellants. 

In  the  simplest  case,  where  only  two  types  of  sites  are  considered 
and  where  the  concentration  of  the  primary  site  is  being  varied, 
the  expression  for  the  detonation  reaction  time  slight  be  described 
in  terms  of  three  different  regions  of  Tr,  viz., 

a.  Primary  Region: 

t1o  <  rr  <  ^20  »  corresponding  to  significant  values  of  N^Q: 
tr  is  simply  given  by  Equation  1-26,  i.e.. 


Ti  „  «  -T. 


Tr  =  TBST - 21”  +  Tlo 

lo 


Tr  =-ri(N^0” '  V 


(1-34) 

(1-35, 


b.  Transition  region: 

Ti0  «  Tr  k  t2o,  corresponding  to  small  values  of  N^o : 
Under  this  set  of  conditions,  Equation  1-28  is  approximated  by 

Tr  =  7®"  [  1  -  hNlo  -  L2N2o  ]  (I- 

Since  L2^2b*  equation  can  be  written  as 

Tr  =“21^7  [1_L2N2o]  CI“ 

c.  Secondary  region: 

Tlo  <  T2o  <  Tr»  corresponding  to  negligible  values  of  Ni0 : 
Again,  Tr  is  simply  given  by  Equation  1-26,  i.e., 


Tr  =1M7’ 


(1-36) 


(1-37) 


Tr  =- 


1 

L2 

1**20 

"  2B  2° 

1  I 

1  T 

a~S W  | 

N  “  “9 

.  w2o  z 

(1-38) 


(1-39) 


CORRELATION  OF  EXPERIMENTAL  CRITICAL  DIAMETER  DATA 

In  Reference  1-7  it  was  shown  by  both  theoretical  and  experimental 
arguments  that  the  critical  diameter,  do,  of  RBX  adulterated  solid 
composite  propellant2  should  be  directly  proportional  to  the  deto¬ 
nation  reaction  tiros  at  the  critical  detonation  velocity,  l.e.. 


dc  *  K(Vc 


(I-W 


‘Convent  tonal  propellant  baaed  upon  AF,  Al ,  and  rubber  binder  in 
which  varying  amounts  of  AS*  were  replaced  by  SDK  (15  u  part  idee). 


where  K  is  a  proportionality  factor  which  remains  relatively  con¬ 
stant  over  an  kua  range  ol  u  to  J.U7« , 

Also  described  in  Reference  1-7  as  well  as  in  Reference  1-9  are 
experimental  critical  diameter  data  for  RDX  adulterated  propellant 
in  which  dc  was  varied  from  ~  2  in.  to  ~  70  in.  as  the  KDX  coi?tent 
was  varied  from  ~  10%  to  0%.  The  reported  critical  diameter  data 
are  shown  in  Table  I- I. 

It  is  interesting  to  examine  these  data  in  light  of  Equation  1-41 
and  the  expressions  for  Tr  which  were  derived  in  the  preceding 
sections. 


Case  of  a  Single  Initiating  Site 

If  the  RDX  particles  are  the  sole  source  of  initiating  sites  then 
the  expression  for  the  reaction  time  becomes  that  given  by  Equa¬ 
tion  1-27.  Con  lining  this  expression  with  Equation  1-40  then 
yields  an  expression  for  dc  of  the  form 


(1-41) 


As  noted  earlier,  Equation  1-41  is  identical  to  that  derived  iri 
Reference  1-7.  A  best  experimental  fit  to  the  data  was  reported 
when  K1  =  15.3  and  K2  =  -30.9  (dc  in  in.).  The  best  fit  of  the 
data  in  Table  I -I  with  Equation  1-41  is  shown  in  figure  1-1.  It 
is  readily  seen  that  in  the  range  of  fRDX  >  0-02,  Equation  1-41 
fits  the  data  quite  well;  however,  at  lower  values  of  fRDX*  the 
theoretical  curve,  (based  upon  the  best  values  of  and  K2)  pre¬ 
dicts  values  of  dc  which  are  too  large.  This  fact  combined  with 
the  fact  that  dc  approaches  infinity  as  fRDX  approaches  zero  led 
0.  R.  Irwin  at  Aerojet  (Reference  1-7)  to  suggest  that, there  were 
additional  initiating  sites  inherent  to  the  propellant  which  act 
as  an  effective  constant  weight  fraction  of  RDX  (see  Equation  1-29) 

In  any  case  it  would  appear  from  Figure  T-l  that  all  the  expei i- 
mental  dc  data  can  not  be  correlated  by  the  Tr  expression  for  a 
single  type  of  initiating  site. 


3This  work  was  carried  out  by  Aerojet-General  Corporation  as 
part  of  the  Air  Force  Solid  Propellant  Hazards  Program  (Pro¬ 
ject  SOPHY) . 
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Table  I-I.  Experimental  Critical  Diameter  Data  for 
RDX -Adult era ted  Solid  Composite.  Propellant 
(References  1-7  and  1-9). 


Mass  Fraction  of 
RDX 

CfRDX) 

Critical  Diameter 
<dc) 

(in.) 

(£rik)1/3 

0.10 

2.12 

0.4642 

0.092 

2.66  /■; 

0  *4514 

0.071 

5.25 

0 .4141 

0.0475 

LI. 25 

0.3621 

0.021 

23.5 

0.2759 

0.00375 

48.0 

0.1554 

0 

/60  (No  Go)\* 

172  (Go)  / 

0 

♦Results  of  one  test  at  each  diameter. 


\\ 


u 


CRITICAL  diameter 


Case  of  x¥o  Types  or  lnitiat ing  bites,  witn 
Negligible  Initiation  times  ” 


In  this  case,  the  one- dimensional  lattice  model  expression  for  'tj. 
(i.e.,  Equation  1-31)  and  Equation  1-40  yields 


K' 


d  = - 

c  (f 


RDX 


1)l/3  +  c,  ~  K' 2  (fRDX) 


1/3 


(1-42) 


The  best  fit  to  the  experimental  data  of  Table  I-I  is  found  for 
the  parameter  values  K'i  =  24.55,  K' 2  =  -50.54,  and  c’  =  0.19. 
The  correlation  of  Equation  1-42  with  the  data  is  shown  in  Fig¬ 
ure  1-2.  Also  shown  is  Irwin’s  correlation  curve  (References 
1-7  and  1-9)  which  corresponds  to  the  use  of  Equation  1-29  for 

Tr ,  i .  e.  , 


(f  +  c)1/3 
URDX  ; 


-  Xr 


(1-43) 


It  is  readily  seen  that  both  expressions  for  dc  correlate  the  data 
very  well  for  fRDX  51  0.00375;  however,  the  results  of  the  two 

experimental  tests  at  zero  percent  RDX  appear  to  favor  Equation 


1-43. 


If  one  assumes  that  the  experimental  dc  data  reflect'  the  existence 
of  inherent  RDX-like  initiating  sites,  then  one  might  conclude 
that  Irwin’s  expression  for  Tr  (Equation  1-29)  may  be  more  appro¬ 
priate  than  the  corresponding  one -dimensional  lattice  model  ex¬ 
pression  (i.e..  Equation  1-33).  It  is  noted  here  that  the  pri¬ 
mary  discrepancy  between  the  two  expressions  lies  in  the  differ¬ 
ence  between  the  factors  (fRDX  4  c)"i/3  and  (fRnx1'3  4  c' )  "S 
which  in  essence  describe  the  average  distance  between  centers 
of  all  initiating  sites.  The  factor  (fRQX^/3  ♦  c')-!  results  from 
averaging  the  one -dimensional  concentration  expressed  as  number 
of  sites/unit  length,  i.e., 

S'*  (zfjV3)"1  (1-44) 

or 

S  oc  (fRcncl/3  +  c*)  1  (1-45) 


4 

i 


! 


1 

1 


1 
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However,  in  three  dimensions  the  average  distance  between  sites 
should  be  based  upon  an  average  of  the  number  of  siter/unit 
volume.  This  would  probably  lead  to 

S  oc  (  2  fp  -1/3  (1-46) 

or 

S  <*  (f^  +  c)"1/3  (  I -47) 

Thus  if  the  assumption  of  inherent  RDX-like  initiating  sites  ia 
valid  then  the  better  fit  of  Irvin's  expression  to  the  experi¬ 
mental  ac  data  would  suggest  that  a  three-dimensional  lattice  model 
may  be  required  for  a  more  quantitative  description  of  detonation 
in  heterogeneous  explosives.  Unfortunately,  the  three-dimensional 
approach  to  the  problem  brings  up  additional  mathematical  orob- 
lems ,  and  is  not  available  at  this  time. 

Alternatively,  the  fact  that  the  apparent  dc  for  £rbx  =  0  falls 
below  the  value  that  is  predicted  by  Equation  1-42 ,  while  the  dc 
data  for  fRiBf  ^  0,02  follow  Equation  1-41  quite  well  (i.e.,  single 
type  of  initiating  site)  might  indicate  that  the  inherent  initi¬ 
ating  sites  are  active  only  for  very  small  values  of  fi*BX  (i.e., 
only  for  large  tr) .  This  brings  us  to  the  possible  situation 
described  earlier  in  this  section  which  is  applied  to  the  dc  data 
in  the  following  paragraphs. 


C tse  Where  the  Number  of  Types  of  Initiating  Sites 
Varies"  with  Tr 

If  one  assumes  that  the  dc  data  of  Table  I -I  reflect  the  exisr 
tence  of  inherent  initiating  sites  which  become  active  only  when 
tr  exceeds  a  certain  value,  then  in  accordance  with  the  approxi¬ 
mation  treatment  shown  earlier  in  this  section  the  data  of 
Table  I--X  should  be  fit  by  three  different  bun  related  expressions 
for  dc  based  upon  three  different  regions  of  T„. 

a.  Primary  region: 

( to)  ppx  -tr  <  Ct0>x 

In  the  primary  region,  only  the  RDX  particles  are  active  initia¬ 
ting  sites,  hence  Equations  1-35  and  1-40  obviously  yield  air  ex¬ 
pression  for  &<;  which  is  the  same  as  Equation  1-41.  Prom  Fig¬ 
ure  .1-1,  it  can  be  estimated  that  the  primary  region  extends  over 
the  range  2.12  =$  dc  <  23.5  in. 
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It  is  now  useful  to  express  the  constants  of  Equation  1-41  in 
terms  of  more  fundamental  quantities.  Noting  that 


N; 


6  P 


it  P.  Li* 


1/3 


«PV3  ■  al£il/3 


( 1*48) 


where  P  is  the  propellant  density.  Equation  1-41  can  now  be 
expressed  as 


(V, 


R 


2R 


c 


a,  f  i 


17J  "Li 


{ 1-49) 


where  the  subscript  1  now  represents  RDX.  Comparison  of  the  con¬ 
stants  of  Equation  1-49  with  those  of  Equation  1-41  yields 


Kj  =  K/2a,Bc  =15.3 
=  KL1/2Bc  =  -30.9 


(1-50) 
( 1-51) 


These  parameters  can  now  be  used  in  the  expression  for  d^  in  the 
transition  region. 

b.  Transition  region: 

'^o'rdX  «  T.-  *  <Vx 

Combining  Equations 1-3 7 ,  1-40,  and  1-48,  yields 


K 


c' 


2a, D  f, 

1  cl 


KL, 


T75  *  2B, 


**20*2 


2Bcalfl 


T7T 


(1-52) 


where  the  subscript  2  refers  to  the  inherent  initiation  sites. 

Utilizing  Equations  1-50  and  1-51,  the  expression  for  d*.  can  be 
expressed  as 


<<U  = 


K 


1  „  V* 


L.  1/  A. 

Ur  -*4  ~-~r/y 

-  n  *i 

where  K3  *  (dimensionless). 


( 1-53) 


Assuming  tha^  dc  =  43  in.  represents  the  transition  region  yields 
a  vain**  Kj  0.1994.  The  resulting  «smpii- ical  reiarionsnip  tor 
dc  in  the,  transition  region  is 


(d  )  =  - -  ~  30.9  (inches)  (1-54) 

transition  ^RDX'’^ 


A  comparison  of  Equation  1-34  with  the  experimental  data  is  shown 
in  I  igure  1-3;  also  shown  is  a  replot  of  Equation  1-41  for  the 
pi imar y  reg ion . 

c.  Secondary  region: 


(  i;  <  (t  )  c  t 

o'RDX  v  o'x  r 

Application  of  Equations  1-39,  1-40, 
yields 


and  1-48  to  the  secondary 


I 


<d').  =  ^  [  K3 

Since  (dc)8  depends  only  upon  the  inherent  ignition  siter,  it  must 
be  a  constant  of  the  propellant  system.  This  is  shown  in  Fig¬ 
ure  1-3  where  (dc)s  is  estimated  from  the  experimental  data  to 
be  70  in.  The  actual  variation  of  dc  with  fRDJf  is  then  assumed  to 
follow  the  dashed  curve  of  Figure  1-3,  which  is  merely  an  esti¬ 
mated  interpolation  of  how  the  three  different  detonation  regions 
transform  into  each  otnei . 

Obviously,  the  approach  taken  above  in  correlating  the  experimen¬ 
tal  data  does  not  offer  the  predictive  characteristics  of  the 
cases  previously  described;  however,  it  does  serve  to  offer  a 
possible  alternate  route  to  explaining  the  detonative  behavior 
of  RDX  adulterated  propellants. 

It  is  interesting  to  note  that  the  above  type  of  data  correlation 
is  not  completely  devoid  of  predictive  capabilities,  since 
choosing  the  constant  (dc)  at  70  in.  sets  the  value  of  a]t,2  at 
1.265  in.  From  Reference  1-7,  can  be  estimated  as  2.07  x  103 
which  yields  a  value  of  1*2  -  6,11  x  10-^  in.  Thus,  assuming  that 
the  ab<-ve  analysis  is  corr  ct,  the  inherent  hot-spot  diameter 
is  ~  16  p.  Also,  since  N2o  =  K3/L2,  the  resulting  distance  be¬ 
tween  hot-spot  centers  (i.e.,  I/N20)  becomes  ~  80  p„ 


-Ki 


(1-55) 
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Figure  1-3.  Correlation  of  Critical  Diameter  Data 

of  Sites  Varying  with  Tr. 


The  grain  burning  distanc  i  between  hot -spots  then  becomes  S  -  Ly 
=  64  p,.  It  is  interesting  to  note  that  in  the  propellant  under 
consideration,  the  AF  uunLeuL  was  composed  or  a  Dimodat  distribu¬ 
tion  in  which  ~  60p  radius  particles  (peak  size)  comprised  ~  50% 
of  the  propellant  mass,  and~  6p  radius  particles  comprised  ~  18% 
of  the  propellant  mass.  It  is  suggestive  from  the  values  found 
for  L2  and  N2o  that  the  inherent  initiation  sites  might  be  the 
small  AP  particles,  which,  along  with  aluminum,  reside  in  the 
interstices  between  the  large  AP  particles.  The  distance  between 
the  interstices  of  the  large  AP  particles  can  be  shown  from  sim¬ 
ple  geometric  considerations  to  be  ~  70p.  It  should  be  empha¬ 
sized  that  at  present  this  suggestion  must  be  considered  highly 
speculative. 


SUMMARY  AND  CONCLUSIONS 

A  lattice  model  of  detonation  of  heterogeneous  explosives  has 
been  described  based  upon  the  concepts  of  both  hot-spot  ignition 
theory  and  grain  burning  theory. 

Through  a  one-dimensional  treatment  of  the  model  it  has  been  pos¬ 
sible  to  derive  the  classical  Eyring  grain  burning  equation  *nd 
to  define  the  assumptions  underlying  the  use  of  the  Eyring  equa¬ 
tion,  viz:  (1)  that  hot-spots  have  no  finite  dimension  or  ig¬ 
nition  time;  and  (2)  the  average  concentration  of  hot-spots  is 
equal  to  the  concentration  of  explosive  granules  (i.e.,  one  hot¬ 
spot  per  granule).  In  addition,  it  has  been  possible  to  relax 
these  assumptions  and  to  extend  the  one -dimensional  treatment  to 
other  explosive  situations,  i.e.,  (1)  where  the  time  to  ignition 
of  grain  burning  is  nonnegligible  compared  to  the  overall  deton¬ 
ation  time,  and  (2)  where  there  may  be  several  types  of  ignition 
sites,  each  type  having  a  discretely  different  ignition  time. 

Such  e.xplosive  situations  appear  to  have  application  to  the  deton- 
ability  of  RUX  adulterated  solid  composite  propellant. 

For  the  case  of  RDX  adulterated  propellant  in  which  the  RDX  par^ 
tides  are  the  sole  source  of  hot-spots,  the  one -dimensional 
treatment  of  the  lattice  model  leads  to  a  detonation  reaction  time 
expression  which  is  identical  to  an  expression  previously  pro¬ 
posed  by  Dr.  Irwin  at  Aerojet.  The  proposed  expression  for  has 
already  been  shown  to  be  very  satisfactory  in  describing  the  vari¬ 
ation  of  critical  diameter  of  composite  propellant  containing 
varying  amounts  of  RDX  in  the  range  of  2  to  10%. 
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For  the  csco  of  props  .lant  co?it“ini.!\g  1  <•«.«*  t-h^n  9%  RDX  the  c.ri f i cal 
diameter  data  suggest  that  RDX  cannot  be  a  so'ife  source  of  hot¬ 
spots.  Extending  the  lattice  model  to  this  case  results  in  ex¬ 
pressions  for  rr  which  differ  somewhat  from  an  expression  which 
was  proposed  and  successfully  applied  by  Dr.  Irwin  for  the  same 
explosive  situation.  However,  the  differences  are  relatively 
slight  and  may  simply  invol'  e  the  use  of  a  one -dimensional  approach 
rather  than  a  three -dimensional  approach  in  describing  the  aver¬ 
age  concentration  of  hot-spots. 

In  any  case,  it  would  appear  that  the  present  lattice  model  of 
detonation  serves  as  a  fundamental  basis  for  the  description  of 
propellant  detonation  which  had  been  initially  developed  at  Aerojet. 

It  would  be  of  future  interest  to  undertake  a  more  rigorous  mathe¬ 
matical  treatment  of  the  lattice  model  in  which  certain  constraints 
(viz.,  those  involved  in  the  approximation  of  Equation  1-19)  were 
removed  from  the  one -dimensional  treatment.  Also,  the  treatment 
of  the  lattice  model  should  be  extended  to  three  dimensions. 

Finally  the  lattice  model  should  be  extended  to  include  the  effects 
of  lateral  quenching  waves'  (i.e.,  x*arefaction  waves)  on  the  over¬ 
all  detonation  reaction  time.  Such  an  extension  could  lead  to  a 
description  of  nonideal  detonation  phenomenon  such  as  fading  deto¬ 
nation  and  charge  diameter  effects. 
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APPENDIX  II 


COMMENTS  ON  THE  RUBY  CODE 

R.  F.  Chaiken ,  Technical  Consultant 
Aero jet -General  Corporation, 
Downey,  California 


INTRODUCTION 

The  RUBY  Code  (Reference  II-l)  is  a  FORTRAN  computer  program 
designed  to  calculate  the  ideal  detonation  properties  of  high 
explosives  utilizing  Chapman -Jouguet  (C-J)  theory.  The  program 
is  based  on  the  assumption  that  the  gaseous  products  obey  the 
Becker -Kistiakowsky -Wilson  (BKW)  equation  of  state  (References 
II-2  and  II-3) ,  which  can  be  written  as 


JKL  - 
nRT 


where 


X  exp  (<9X) 
<5-i  ^ini 


(II-l) 


V  (T  +8)<*  n± 

Here,  a,  8,  x,  9*  ki  ere  constants,  and  hi  ia  the  mole  fraction 
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gaseous  product. 


In  RUBY,  this  P,  V,  T  relation  is  used  to  express  the  one-dimen¬ 
sion  detonation  conservation!  equations  and  C-J  hypothesis  (i.e., 
D  *  Cj  +  Wj)  in  terms  of  P,  V,  T  variables,  and  also  to  describe 
the  fugacity  of  the  gaseous  reaction  product  2  cs  a  function  of 
pressure  and  temperature.  In  addition,  RtL?  employs  an  equation 
of  state  c°  the  form 
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to  handle  the  possible  existence  of  one  or  two  solid  detonation 
products. 


The  T-J  t,tate  is  assumed  to  he  «t  t n Krnmfl yn. amic  equilibrium  (isini- 
mum  i.ree  energy)  and  is  calculated  by  a  method  of  steepest  descent 
which  is  described  by  White,  et  al  (Reference  II-4;. 

The  main  use  of  RUBY  to  date  has  been  in  the  calculation  of  ideal 
detonation  properties  of  CHNO  explosives  (References  II-2,  II-3, 
II-5,  11-6)^.  The  general  approach  has  been  to  determine  the  BKW 
constants  'i.e,,  a,  6,  x,  9,  k^)  which  allow  the  best  fit  to  the 
experimental  D(Pg)  and  Pj(PQ)  data  for  one  or  two  specific  ex¬ 
plosives  (e.g.,  RDX,  PETN) ,  and  after  selecting  these  constants, 
to  apply  RUBY  to  other  explosives. 

)  for  CHNO  explo- 
owever,  there  is 

apparently  no  single  set  of  BKW  parameters  which  yields  good  agree 
ment  with  all  20  of  the  explosives  considered.  For  example. 

Table  4  in  Reference  II-5  illustrates  that  detonation  velocities 
and  pressures  may  be  in  error  by  as  much  as  10%  and  15%,  respec¬ 
tively.  Also,  while  there  is  a  lack  of  reliable  experimental  data 
on  detonation  temperatures,  it  would  appear  that  RUBY  values  of 
Tj  may  be  up  to  40%  too  low. 

The  purpose,  of  this  report  is  to  examine  the  results  of  the  RUBY 
program  and  to  determine  its  usefulness  in  calculating  the  ideal 
detonation  properties  of  conventional  solid  propellants  containing 
ammonium  perchlorate,  aluminum,  and  oxygen-deficient  rubber-type 
binders.  Of  particular  interest  is  the  use  of  RUBY  to  determine 
the  effect  of  incomplete  chemical  reaction  on  the  ideal  detona¬ 
tion  properties  (D,  Pj,  Tj).  Toward  this  end,  RUBY  calculations 
have  been  carried  out  for  ammonium  perchlorate  CAP),  alone  and  in 
combination  with  typical  propellant  ingredients,  and  for  RDX/ 
aluminum  mixtures .  Computer- input  techniques  were  developed  to 
allow  various  amounts  of  aluminum  to  remain  unreacted.  These  cal¬ 
culations  brought  out  certain  apparent  internal  inconsistencies 
which  suggest  that  RUBY  should  not  be  used  to  predict  the  detona¬ 
tion  properties  of  aluminized  explosives. 


Generally,  the  RUBY  calculated  D(PQ)  and  Pj(P 
sives  are  in  fair  agreement  with  experiment;  I 


RUBY  CALCULATIONS 


Ammonium  Perchlorate 

The  detonation  properties  of  ammonium  perchlorate  (AP)  have  been 
calculated  by  RUBY  over  the  range  of  loading  densities  0.8  ^  P0 
1.5  gm/cc.  Various  sets  of  BKW  constants  have  been  employed  to 


^References  11-2  and  II -3  do  not  refer  to  RUBY  calculations  per 
se,  but  to  similar  calculations  using  tba  BKW  equation  of  state. 


attempt  to  match  the  reported  experimental  data  (References  II-7 
and  II-8) ,  as  well  as  to  ascertain  the  sensitivity  of  the  com¬ 
puted  detonation  properties  to  changes  in  the  BKW  constants. 

The  results  of  these  calculations  are  shown  in  Figures  11-1  through 
11-3  where  ideal  detonation  velocity  (D),  Chapman -Jouguet  pres¬ 
sure  (Pj),  and  Chapman -Jouguet  temperature.  (Tj)  are  plotted  against 
Pq.  The  notes  for  Figures  II-l  through  II-3  describe  the  con¬ 
ditions  for  obtaining  curves  A  through  F.  These  results  show  that 
by  suitable  adjustment  of  the  BKW  constants  (a,  x,  B,  and  in 
particular,  kj  for  the  principal  gas  products,  e.g. ,  HC1)  almost 
any  linear  D(P0)  curve  can  be  obtained. 

However,  the  constants  that  have  been  derived  for  best  fit  with 
CHNO  explosives  (e.g.,  curve  A  and  B) ,  although  yielding  reason¬ 
able  agreement  with  the  experimental  Pj  value  at  Po  =  1.0  gm/cc , 
do  not  yield  good  agreement  with  the  experimental  D  value.  Even 
curve  E,  which  is  presumably  the  result  of  SRI's  attempt  (Refer¬ 
ence  II-8)  to  optimize  the  BKW  constants  for  AP,  falls  short  of 
being  in  good  agreement  with  experimental  detonation  velocities. 

With  regard  to  the  RUBY-calculated  Chapman -Jouguet  temperatures , 
it  is  readily  seen  that  an  increase  in  PD  results  in  a  decrease 
in  Tj.  This  Tj(P0)  relationship  appears  to  be  common  to  all 
RUBY  calculations,  including  those  for  CHNO  explosives  (Refer¬ 
ences  II -2,  II-3,  II-5,  and  II -6).  This  undoubtedly  arises  from 
the  fact  that  the  BKW  equation  of  state  considers  only  a  repul¬ 
sive  potential  between  the  detonation  products.  On  the  other  hand, 
curve  F,  which  corresponds  to  calculations  with  Cook's  covolume 
equation  of  state  (Reference  II-9) ,  shows  Tj(P0)  to  be  an  in¬ 
creasing  function  of  loading  density.  Cook's  covolutne  depends 
only  on  volume  and  hence  does  nol  consider  the  potential  energy 
arising  from  intertnolecular  interaction.  Unfortunately,  the 
present  lack  of  reliable  experimental  Tj  data  precludes  a  clear- 
cut  answer  as  to  whether  Tj  should  be  an  increasing  or  decreasing 
function  of  p0. 

In  the  case  of  high-density  CHNO  explosives,  it  has  been  argued 
(Reference  II-2)  that  since  the  C-J  density  (Pj)  is  generally 
greater  than  the  explosive  crystal  density,  the  distances  between 
atomic  and  molecular  species  In  the  C-J  plane  are  so  small  so  that 
the  interactions  between  the  species  are  primarily  repulsive 
(hence,  the  BKW -type  of  equation  of  state). 

In  this  case,  even  though  the  total  change  in  specific  internal 
energy  of  detonation  increases  with  PQf  the  net  result  could  be 
to  increase  the  potential  energy  of  the  C-J  system  at  the  expense 
of  the  kinetic  energy. 


tun  veLoarr 


I  I  I  I  I 

Curve  At  Covolum*  set  ueed  by  t>Ri  In  calculating  Dfp0)  for  AP 

(Privets  communication  from  L.B.  r>M l«y  (January  4 966)). 

Curve  C:  Cavolu.vte  sat  used  by  NOTS  in  calculating  propellent 

detonation  proportion  (Priyata  comaiunlcction  front  J  .  Dleroff 
(Mai  Lb  1964)}. 
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Curve  C:  Cove  fume  Mf  by  NOTS  in  cakultHnf  propellent 

detonation  properties  (Private  ctMmnwnicatioo  from  J.  Dior  off 

(March  1964))*' 

Curve  D  i  Arbitrary  dec  root#  in  eevoluiMt  of  MCI  and  Clj  (raw  those 
of  Ciau  C, 

Curvo  E:  Sfll  RUBY  calculations,  from  Reference  )|  .g, 

Curvo  P;  A&C  calcvlotiMi,  from  Refers  co  II  -7#  using  Cook's 

covolwmo  equation  of  state  fitted  to  the  experimental  Ideal 
detonation  velocity  at  p  =s  1.0  gn/Cc,* 
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Figure  I I -3.  Calculated  Detonation  Temperatures 
of  Ammonium  Perch lor at a. 


NOTES  FOR  FIGURES  I I -I  THROUGH  I I -3 


Curves  A  and  B 

RUBY  calculations  to  determine  the  effect  of  a  change  in 
BKW  constants,  using  the  same  thermodynamics  data: 

Curve  A:  a  =  0.5;  3  *  0.09;  x  «  11.85;  6  =400 

Curve  B:  a  =  0.5;  3  =  0.16;  x  =  10.91;  0  =  400 

The  BKW  set  for  curve  A  corresponds  to  the  set  used  by 
Cowan  and  Fickett  (Reference  II-2)  for  curve-fitting 
65/35  RDX/TNT  D(Pq)  data.  The  BKW  set  for  curve  B  cor¬ 
responds  to  the  set  used  by  Mader  (Reference  11-3)  for 
curve  fitting  RDX  D(P0)  data. 


Curves  A,  C,  and  D 

RUBY  calculations  to  determine  the  effect  of  a  change 
in  covoliune  constants  using  the  same  BKW  constants. 
Values  of  kf  for  curves  A,  C,  and  D: 


Curve  C 

Curve  D 

Curve  A 

h20 

360 

360 

250 

h2 

180 

180 

180 

N2 

380 

380 

380 

nh3 

476 

476 

476 

n20 

670 

670 

670 

NO 

336 

386 

386 

°2 

350 

350 

350 

no2 

670 

670 

600 

HC1 

1588 

794 

643 

ci2 

1157 

578 

532 

However,  for  ■mmon i utn  perchlorate  at  the  loading  densities  con¬ 
sidered  here..  P  i  i  n  less  than  the  crystal  density  M  Qh  rrm/rr) 
it  is  difficult"1  therefore  to  understand  why  the  above  argument 
should  still  be  valid  under  these  conditions.  It  is  believed 
that  the  RUBY  -calculated  decrease  in  Tj  as  p0  increases  is  un¬ 
realistic  for  AP.  On  the  other  hand,  it  should  be  stated  that 
the  Tj  values  obtained  for  AP  by  use  of  the  Cook-covolume  equa¬ 
tion  of  state,  which  are  600  to  700°K  higher  than  the  AP  adia¬ 
batic  flame  temperature  (1  atm),  might  also  be  unreasonable. 
However,  the  increasing  Tj(P0)  function  that  is  obtained  by  such 
a  covolume  equation  of  state,  is,  at  least  to  this  writer,  intu¬ 
itively  plausible. 


Solid  Composite  Propellant 


RUBY  calculations  have  been  carried  out  for  a  solid  composite  pro¬ 
pellant  composition  containing  AP ,  aluminum,  and  PBAN  binder 
( polybutadiene-acrylonitrile  copolymer).  The  effect  of  replacing 
part  of  the  AP  with  RDX  on  the  detonation  properties  has  been 
calculated,  as  well  as  the  effect  of  nonreaction  of  the  aluminum. 
This  latter  effect  corresponds  to  a  current  theory  that  the  alumi¬ 
num  oxidation  reaction  proceeds  too  slowly  to  occur  within  the 
detonation  reaction  zone  (Reference  11-10). 

Prevention  of  aluminum  oxidation  in  the  RUBY  calculation  is  readily 
accomplished  by  replacing  all  or  part  of  the  aluminum  in  the  ex¬ 
plosive  composition  by  the  fictitious  metal  AlX.  The  AlX  has  all 
the  thermodymudc  properties  of  Al,  but  no  oxidation  products 
(e.g.,  AIX2O3,  A1X02,  AIXCI3). 

The  results  of  the  RUBY  calculations  with  the  BKW  constants  and 
thermal  data  corresponding  to  curve  G  in  Figures  II-l  through 
1,1-3  ar  shown  in  Table  II-I.  The  data  fo~  propellant  A  (normal 
case,  w  th  100%  Al  reaction)  are  comparable  to  data  obtained  for 
high  explosives  with  similar  heats  of  explosion.  A  comparison 
of  propellant  G  (9.2%  RDX  adulterated  propellant)  with  propel¬ 
lant  A  shows  that  the  detonation  parameters  (E,  Pj,  Tj,  Ej-Ec) 
all  increase  when  the  more  energetic  RDX  (  AQ  »  1300  Kcal/gm)  re¬ 
places  AP  ( aQ  -  300  Kcal/gm).  These  results  are  to  be  expected. 

The  effect  of  replacing  aluminum  (reactive)  with  AlX  (unreactive) 
is  that  Pj,  Tj,  Ej-Ec  decrease  as  might  be  expected,  but  there  is 
a  4%  increase  in  the  detonation  velocity.  This  is  evident  by 
comparing  propellant  A  with  propellant  B. 
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Table  II-I.  Detonation  Properties  of  AP  Propellant  as 

Calculated  by  SU5Y. 

Propellant  A:  AP/Ai/PBAA  =  0.69/0.15/0.16  (1.00%  A1  reaction) 

TropelLant  B:  AP/A1/P8A4  =  0.69/0.15/0.16  (5%  Al.  reaction)* 

Propellant  C;  AP/A1/PBAA/RDX  *  0.598/0.15  Q. 16/0. 092.  (100%  Al 

reaction) 


DETONATION  PROPERTIES 

A 

P  =1.73 

0 

Propellant 

B 

P  =■  1.73 
o 

c 

p0  =  1.715 

D,  mm/psec 

7.14 

7.40 

7.25 

Pj,  Kbar 

226 

206 

232 

Tj,  PK 

3198 

1091 

3217 

Vj ,  cc/mole  of  gaa 

12.39 

14.41 

12.61 

P©/P  j 

0.743 

0.783 

0.743 

Ej  -  Eo,  eal/gm  HE 

401 

308 

415 

BKU  parameter 

5.71 

11.38 

5.86 

_3 

G-J  compoaiticn*  10  molca/gm  HE 

Total  Gases  (♦,**) 

28.85 

27.74 

28.67 

ch4 

5.19 

5.69 

5.76 

CO 

4.25 

nil 

4.70 

GO, 

1.65 

5.40 

1.88 

Cl 

5.73 

nil 

4.91 

CIO 

nil 

5.04 

nil 

H2 

0.19 

nil 

0.15 

h2o 

8.44 

8.29 

6.95 

*2 

2.61 

2.83 

3.44 

NH- 

0.65 

0.20 

0.69 

AlCl 

0.14 

nil 

0.16 

AlCl3 

nil 

0.28 

nil 

Total  Solids  <♦,**) 

2.71 

5.28 

2.69 

AH1) 

— 

5.28 

— 

ai2o3 

2.71 

nil 

2.69 

C  (graphite) 

oil 

nil 

<*)  Unreacted  Al  goes  to  Al(l)  ae  «  detonation  product. 
(♦)  Only  product*  10"^  mole*/^i  H5  are  included. 

(**)  A  daeh  indicate*  the  product  waa  not  prograoeecd. 
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Examination  of  the  gaseous  detonation  products  shows  that  the  total 
amount  of  gas  is  che  same  in  the  two  cases,  and  that  except  foj  an 
increase  in  C02  with  an  accompanying  decrease  in  CO,  there  As  little 
change  in  gas  products  when  aluminum  does  not  react.  From  the 
heats  of  formation  of  AI2O3 ,  CO,  and  COo ,  the  net  enthalpy  loss 
from  the  explosive  system  with  A1X  would  be  ~  350  cal/gm  of  ex¬ 
plosive.  It  vould  be  difficult  to  reconcile  an  increase  in  deto¬ 
nation  velocity  with  this  energy  loss.  Undoubtedly,  the  calculated 
increase  in  D  stems  from  the  value  of  the  BKW  parameter  (X  in 
Equation  II-l) ,  which  does  change  appreciably.  Since  for  CO2 
is  670  vs  390  for  CO  in  these  calculations,  it  is  seen  that  the 
value  of  k  ^ikjm  in  Equation  II-l  will  increase  when  CO  is  con¬ 
verted  to  CO2.  (When  aluminum  does  not  react,  there  is  more  oxy¬ 
gen  available  in  the  explosive  for  reaction  with  carbon.  Thus, 
a  greater  portion  of  the  carbon  is  converted  to  CO2  than  would  be 
the  case  if  aluminum  reacted  to  AI2O3.) 

When  this  increase  is  combined  with  a  decrease  in  C-J  temperature 
(which  results  from  the  A1X  not  reacting) ,  the  BKW  parameter  be¬ 
comes  excessively  large.  It  is  apparent  that  the  accompanying 
increase  in  Vj  (which  results  from  a  lower*  Pj)  is  far  from  enough 
to  prevent  the  excessively  high  value  of  the  BKW  parameter.  The 
apparent  result  is  an  increase  in  the  detonation  velocity,  which 
conflicts  with  the  general  expectation  that  B  is  an  increasing 
function  of  Ej-Eq  for  any  given  explosive.  Similar  results  are 
obtained  for  RUBY  calculations  with  RDX/A1  mixtures. 


Aluminized  Explosives 

Table  II -II  shows  the  results  of  RUBY  calculations  for  RDX  in 
combination  with  Al,  A1X,  and  AI2O3  (80/20  mixtures)  at  the  same 
loading  density  (P^  =  1.94  gm/cc).  For  additional  comparison, 
the  results  of  a  similar  calculation  for  RDX  at  a  slightly  differ¬ 
ent  density  (PQ  =  1.8  gm/cc)  are  also  shown.  These  calculations 
were  carried  out  with  the  same  set  of  thermodynamic  data  and  BKW 
constants.  (The  set  of  input  data  are  thot  1  obtained  frdm  UCRL 
except  for  the  addition  of  input  data  required  for  the  aluminum - 
containing  products:  A10(,  ),  Al20(g),  Al202(g) ,  Al(l),  and 
AI2O3VC).  Presumably,  the  UCRL  data  set  has  been  optimized  for 
RDX.  ) 

While  a  direct  comparison  of  the  calculated  properties  of  explo¬ 
sive  A  (pure  RDX)  with  explosive  B  (20%  aluminum)  is  not  possible 
because  the  loading  densities  are  not  the  same  in  both  cases,  the 
addition  of  reacting  aluminum  greatly  increases  the  detonation 
temperature  (ATj  «  2000°K).  This  is  further  borne  out  by  the 
drastic  drop  in  Tj  when  the  aluminum  is  prevented  from  reacting 
( explosive  C) . 


Tahl  p  TI-II.  Detonatiorii  Properties  of 
as  Calculated  by  RUBY. 

Explosive  A:  RHX 

Expli  * iv*  B:  RDX/A1  -  80/20  (100%  A1  reaction) 
Explosive  C:  REK/A1X  =  80/20  (0%  A1  reaction)* 
Explosive  Ds  RDX/A1203  =  80/20 


RJDX  Explosives 


DETONATION  PHOPERTIES 

A 

Expl< 

» 

waives 

C 

D 

Pq,  gra/cc 

1.80 

1.94 

1.94 

1.94 

D,  ram/^aec 

8  .03 

8.56 

9.11 

8.28 

PJ(  Kbar 

287 

359 

332 

297 

Tj,  °K 

1824 

3882 

2262 

2428 

Vj  P  cc/mole  of  gas 

11.71 

12.28 

11.38 

12.05 

N 

O 

CL 

0.752 

0.747 

0.794 

0.776 

j  Ej  -  Eo,  cal/gm  HE 

472 

464 

421 

410 

BKW  parameter 

9.23 

• 

6.86 

8.71 

8.09 

C-J  composition,  10”^  molc^/gm  HE 

Total  Gases  (*,**) 

33.80 

25.32 

27.06 

27.21 

ch4 

0.38 

3.33 

0.04 

0.36 

CO 

0.10 

4.47 

0.17 

0.57 

C02 

7.09 

2.48 

5.45 

5.73 

H2 

nil 

0.04 

nil 

nil 

h2° 

12.72 

2.85 

10.53 

9.57 

«2 

13.50 

10.38 

10.74 

10.64 

nh3 

0.01 

0.83 

0.13 

0.34 

NO 

nil 

0.05 

nil 

nil 

A120 

— 

0,91 

-- 

nil 

Total  Solid*  (♦,**) 

5.92 

3.31  ' 

12,55 

6.09 

aku 

-- 

— 

7.41 

Al203  I 

-- 

2.80 

--  • 

1.96 

G  (graphite) 

_ 

0.51 

5.14 

4.13 

(*)  Li  art  acted  A1  got*  to  Aid)  ee  a  detonation  product, 
<  ♦)  Only  product*  >  IG_S  molea/ga  HE  are  included. 

<**)  A  deeh  indicate*  the  product  woe  not  programed. 
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As  in  the  aluminized  prop®!!®™*-  rase,  it  is  readily  seen  that  non- 
reaction  of  aluminum  (explosive  G)  causes  an  increase  in  D,  while 
Pj ,  Tj,  and  Ej  all  decrease.  Again,  the  relatively  lar^e  increase 
in  the  BKW  parameter  with  A1X  suggests  that  the  conversion  of  GO 
to  GO2  may  be  responsible  for  this  effect. 

However,  it  appears  that  replacing  aluminum  by  Al^Oj  (explosive  D) 
does  not  show  a  similar  effect,  even  though  in  this  case  AI2O3 
can  be  considered  as  an  inert  ingredient  in  the  same  manner  as 
A1X.  A  comparison  of  the  detonation  products  of  explosives  G 
and  D  does  not  appear  to  offer  any  significant  clues  to  the  dif¬ 
ferent  effects  of  the  two  inert  additives  on  RDX.  However,  the 
effective  heats  of  formation  used  for  the  two  RDX  mixtures  were 
somewhat  different:  76.8  cal/gm  for  explosive  C  vs  686  cal/gm 
for  explosive  D.  It  is  suggested  that  under  these  circumstances 
a  direct  comparison  of  the  two  cases  may  not  be  too  meaningful 
without  experimental  data  for  both  RDX/A1  and  RDX/AI2O3  mixtures. 

There  are  abundant  experimental  data  on  aluminized  high  explosives 
which  show  that  aluminum  lowers  the  detonation  velocity  of  the 
pure  explosives  (Reference  II-9) .  It  has  been  proposed  that  the 
aluminum  either  behaves  as  an  inert  diluent  or  that  it  reacts  to 
Al20(g)  with  an  overall  endothermic  effect.  The  results  of  the 
RUBY  code  appear  to  be  inconsistent  with  the  experimental  results 
as  well  as  with  either  of  these  explanations. 


DISCUSSION  AND  CONCLUSIONS 

From  the  RUBY  calculations  carried  out  for  CHNO  explosives  and  for 
AP  propellants,  it  is  apparent  that  with  sufficient  adjustment  of 
the  many  available  parameters  in  the  BKW  equation  of  state,  suit¬ 
able  D( PQ)  and  P(PG)  data  can  be  calculated  for  any  given  explo¬ 
sive.  However,  it  is  likewise  apparent  that  the  extension  of  the 
RUBY  calculations  to  other  explosives  with  selected  BKW  constants 
may  lead  to  highly  questionable  results,  particularly  when  those 
explosives  involve  new  gaseous  products.  Although  this  point  had 
already  been  made  clear  by  Cowan  and  Fickett  in  their  original 
paper  (Reference  II-2) ,  it  has  apparently  not  been  emphasized  by 
subsequent  investigators  who  have  reported  RUBY  calculations. 

In  spite  of  any  agreement  that  one  can  obtain  between  RUBY -calcu¬ 
lated  and  experimental  detonation  velocities  and  pressures,  there 
are  certain  internal  inconsistencies  wt>ich  throw  doubt  on  the  use¬ 
fulness  of  the  BKW  equation  of  state  in  calculating  detonation 
properties.  These  internal  inconsistencies  involve  (1)  a  calcula¬ 
ted  C-J  temperature  which  consistently  decreases  as  P0  increases 
for  any  given  explosive,  and  (2)  the  increase  in  detonation  velocity 


when  tr he  aluisiuuiu  in  alumiuitcu  cAjjtuH ives  uuet>  not  react  to  con¬ 
tribute  energy  to  the  G-J  state.  Both  of  these  effects  appear  to 
be  related  to  the  BKW  equation  of  state  which  considers  the  energy 
of  molecular  interaction  to  be  solely  repulsive.  Under  the  deto¬ 
nation  conditions  which  are  normally  calculated  by  RUBY,  it  is  felt 
that  such  an  equation  of  state  will  overemphasize  the  role  of  pres¬ 
sure  and  underemphasize  the  role  of  temperature  in  determining  the 
product  distribution  and  the  detonation  properties. 

In  this  connection,  it  is  interesting  to  examine  the  actual  value 
of  the  apparent  molar  covolume  (i„e.,  the  excluded  volume)  that 
results  from  the  RUBY  calculations  presented  for  propellants  in 
Table  II-X.  By  comparing  Equation  II-l  with  a  covolume  equation 
of  state,  i.e., 

P(V-b)  =  nRT  ( II-3) 

it  can  be  readily  shown  that  the  apparent  molar  covolume  b/n  can 
be  expressed  as 

b/n  =-B£-  x  exp  flX  (II-4) 

For  propellants  A,  B,  and  C  the  values  of  b/n  are  11.0,  13.5,  and 
11.3  cc/mole,  respectively.  These  values  represent  approximately 
90  to  95%  of  the  calculated  molar  gas  volume  (Vj),  thus  indicating 
an  extremely  compact  G-J  state.  It  seems  unlikely  that  such  a 
state  can  exist  and  still  be  composed  of  recognizably  independent 
molecular  species. 

Although  the  RUBY  calculations  that  have  been  carried  out  to  date 
do  not  exclude  the  possibility  that  the  noted  inconsistencies  can 
be  resolved  by  a  complete  change  in  the  BKW  constants,  they  do  in¬ 
dicate  the  unsatisfactory  nature  of  the  present  code  with  regard 
to  calculating  detonation  properties  of  propellants  and  aluminized 
explosives.  In  view  of  this,  and  the  arguments  presented  against 
an  equation  of  state  that  implies  solely  repulsive  forces,  it  is 
suggested  that  further  work  with  the  RUBY  code  not  be  continued. 

Pending  development  of  a  more  satisfactory  computer  program  (pre¬ 
sumably  incorporating  a  more  realistic  product  equation  of  state) , 
it  is  proposed  that  the  Parlin-Andersen-Miller  procedure  (Refer¬ 
ences  11-11  and  11-12)  used  in  the  SOPHY  I  program  continue  to 
l®  used  to  estimate  the  ideal  detonation  velocity. 
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Figure  III-2,  Side-On  Overpressure ,  Station  6,  Test  CD-96, 
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Figure  1II-4.  Side -Or  Overpressure,  Station  8,  Test  CD-96 
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Figure  III-6.  Face. -Or.  Overpressure ,  Test  CD-96. 
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APPENDIX  IV 


JETTING  Tp, noMEN ON 


Streak-camera  observation;:  of  previous  critical  geometry  tests 
of  hollow  cylindrical  samples  showed  that  an  abrupt  cessation  of 
an  apparently  steady  detonation  occurs  at  some  point  along  the 
length  of  the  sample  (Reference  IV~1).  The  location  of  this 
point  appeared  to  be  determined  by  the  web  thickness  of  the  grain 
when  the  results  obtained  for  only  one  propellant  formulation 
were  analyzed.  Since  it  is  observed  that  a  high-velocity  jet  is 
produced  in  the  core  ol  an  end-initiated  hollow  cylinder,  the 
responsibility  for  the  ar.amalous  behavior  of  hollow  cylindrical 
samples  may  be  ascribed  to  effects  caused  by  the  jetting  pheno¬ 
menon.  Therefore,  a  series  of  tests  was  designed  to  investigate 
the  manner  in  which  jetting  aff^^r®  the  detonation  reaction  of 
hollow  cylinders.  These  tests  include  studies  of  the  effect  of 
core  diameter  (ID),  sample  length,  and  web  thickness  on  the  be¬ 
havior  of  hollow  cylinders. 


EFFECT  OF  CORE  DIAMETER 

Hollow  samples  cast  from  AAB-3139  propellant  having  web  thickness 
equal  to  1.4  in.  were  prepared  in  several  sizes  in  which  only  the 
ID  varied.  Core  diameters  of  0.06,  0.12,  0.25,  0.50,  0.80,  1.5, 
and  3.0  in.  were  chosen.  All  sample  lengths  were  four  times  the 
OD.  The  samples  were,  instrumented  with  four  sets  of  ionization 
probes  in  circuit  with  four  rasteroscillographs  to  obtain  dis¬ 
tance-time  data,  on  the  reactive  shock  wave  and  the  jet.  Three 
sets  of  probes  placed  in  the  propellant  gather  the  data  at-  -radial 
positions  1/4- ,  1/2- ,  and  3/4-web  in.  from  the  outside  of  the 
samples.  The  foux-th  set  of  probes  monitors  the  travel  of  the  jet 
along  the  charge  axis.  Figure  IV-1  shows  the  test  setup. 

A  1-in.  thick  Plexiglas  plate  is  placed  between  the  Composition  B 
booster  and  the  propellant  acceptor.  The  Plexiglas  attenuates 
the  velocity  of  the  shock  wave  entering  the  acceptor.  This  allows 
better  resolution  of  the  initial  behavior  of  the  reactive  shock 
wave  in  the  propellant  by  preventing  its  being  masked  by  the  nor¬ 
mal  attenuation  of  a  highly  overboostered  initiating  wave  in  the 
acceptor,  which  would  occur  if  the  booster  cnarge  were,  placed 
directly  on  the  propellant  sample.  The  Plexiglas  barrier  also 
prevents  the  booster  shock  wave  from  forming  a  jet  in  the  hollow 
acceptor  before  the  propellant  detonation  reaction  produces  its 
jet.  The  1-in.  thick  attenuator  does  not  reduce  the  shock  pres¬ 
sure  below  the  minimum  level  required  to  initia  e  detonation  in 
these  samples. 
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Figure  IV-1.  Jetting  Phenomena  Teat  Setun. 


The  spicific  purpose  of  these  tests  was  to  determine  the  influence 
of  the  ID  on  the  location  of  the  point  along  the  charge  length 
where  detonation  ceases.  it  is  of  interest  to  learn  whether  a 
minimum  core  diameter  exists  below  which  the  jet  has  no  effect  on 
the  propellant  behavior,  and  whether  a  maximum  core  diameter  of 
reasonable  proportions  exists  above  which  the  jet  has  no  effect 
on  the  propel lan l  behavior* 

The  dx stance-vs-time  data  were  reduced  for  plotting  the  detonation 
velocity  vs  distance  along  the  sample  at  each  radial  distance  in¬ 
to  the  sample. 

Figure  IV- 0  is  a  plot  of  the  data  obtained  at  the  l/4~web  depth. 
Only  straight  lines  were  used  to  connect  the  individual  data 
points,  because  to  have  fit  the  data  to  smooth  curves  would  have 
made  it  more  difficult  to  isolate  the  data  of  any  one  sample  sine 
from  the  others.  From  these  data  it  is  evident  that  regardless 
of  the  size  of  the  ID  the  detonation  began  to  fade  after  7  in.  of 
travel  in  the  samples. 

The  results  obtained  at  1/2 -web  depth  are  shown  in  Figure  IV-3 
and  are  consistent  with  the  preceding.  The  dat.'  indicate  that 
from  the  middle  of  the  web  toward  the  outer  surface  of  the  hollow 
samples  the  detonation  wave  front  is  perpendicular  to  the  charge 
axis  until  rapid  fading  begins. 

The  data  obtained  at  the  3/4-web  deptn  (Figure  IV-4)  are  much 
more  erratic  than  those  from  locations  more  removed  from  the  in¬ 
ner  surface.  They  also  show  no  correlation  with  core  size,  but 
do  indicate  that  fading  begins  somewhat  sooner  at  this  position, 
i. e. ,  at  5  to  6  in. 

The  jet-velocity  data  (Figure  IV-5)  show  some  oscillatory  be¬ 
havior.  Moi'e  important,  there  is  evidently  a  correlation  between 
the  size  of  the  perforation  and  the  velocity  of  the  jet.  It  was 
observed  that  higher  jet  velocities  occurred  at  smaller  c  ire  sizes. 
This  would  be  expected  because  the  Mach  interaction  that  produces 
the  jet  would  be  greatest  when  i he  core  size  is  minimum. 

Some  experimental  difficulty  in  probe  placement  with  the  0 . 06 - in . ~ 
ID  samples  prevented  jet -velocity  data  from  being  obtained  beyond 
the  4~in,  distance.  Probes  further  down  the  charge  indicated  low 
velocities  more  typical  of  the  values  obtained  within  the  web, 
and  these  were,  assumed  to  have  been  owing  to  failure  to  reach  the 
small,  core  with  the  probes. 

From  these  tests,  it  is  concluded  that  the  abortion  of  sustained 
detonation  velocity  is  independent  of  core  size  over  the  50-fold 
range  from  0.06-in.  to  3.0-in.  ID. 
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0 ,25-in.  ID 


EFFECT  OF  LENGTH 


To  determine  whether  the  abortion  of  steady-state  detonation  is 
caused  by  an  end-effect  mechanism,  it  is  necessary  to  test  sam¬ 
ples  that  vary  in  length.  In  the  series  of  tests  described  in 
tiie  paragraphs  entitled  Effect  of  Core  Diameter,  samples  having 
lengths  equal  to  four  times  the  OD  were  used.  Since  the  web 

thickness  was  held  constant  the  OD  varied  in  the  same  inci ements 

as  the  ID,  and  the  lengths,  therefore,  varied  over  a  two-fold 
range.  Since  no  correlation  was  found  between  the  location  of 
the  fadeout  point  and  the  ID  it  follows  that  no  correlation  ex¬ 
isted  with  the.  sample  lengths  either.  However,  a  series  of  tests 
planned  explicitly  to  in  estigate  the  length  effect  using  samples 
of  identical  cross-tic  t  ion  < 4.50-:in.  OD  by  1. 50-in.  ID),. 

Four  sample  lengths  were  tested  18  in.  ( 4  x  OD)  22.5  in.  (5 
x  OD)  ,  27  in.  (6  x  OD)  ,  and  J6  in.  ( 1  x  Oi  ).  The  probe  instru¬ 
mentation  was  limited  to  two  Tenths:  one  set  inserted  to  1/4- 

web  thickness,  the  other  io  3,  -web  thickness.  It  xiglas  was 
used  to  attenuate  th»  booster  shock,  wa  2. 

Figure  IV-6  shows  the  da1  a  obtained  from  .he  ul  e  set  of  u»obe;, 
reduced  to  detonation  velocity  vs  distance  along  he  s.mple.  It. 
is  evident  from  these  data  that  the  sudue  t  fading  of  i  ie  detona¬ 
tion  velocity  is  not  caused  by  an  enc  effect.  The  fading  begins 
after  approximately  24  in.,  and  it  is  therefore  n  observed  r 
the  two  smaller'  Lengt  ;.  It  should  be  noted  in  romoai  ing  he  r  > 

suits  for  the  1.5-in,  web  thickness  exu  nd.  thr  dstan  ■  along 

the  sample,  through  •••hich  a  su:  tain  d  dt  r  natioi  veloci  y  iccurs, 
from  7  to  24  in  Thcs  shows  proof  tfiat,  or  i  g  ven  mater  ial,  the 
web  size  is  the  principal  influence  on  the  ’uratica  of  a  sustained 
detonation  velocity,  as  judged  oy  velocity  a  a  obtained  near  tie 
outer  surface  t  ?  ho.Liow  samf  les. 

Figure  IV-7  reveals  '  h.  xt  rat  rdinary  det>  n,:ti<  vert  tty  data  ob¬ 
tained  from  the  inn  i  se  t  of  probes.  The  'eloc  t^  >.s  i^rea;  d  tr 

near  7.5  mm/us<  at  toe  .  4  -  in .  distance,  \gatn,  »  be  vior  >t 
the  detonation  ne..  the  uncr  surface  o  i  ollow  i  ,  Linders  is  etr 
to  differ  sig  lif  <  eu.  ciy  j  rom  ttat  near  the  >  ter  urface.  The 
outer  portion  ol  the  web  behaves  more  t  "'mai  i  y  unt  i  l  t  re  p.  int  is 
readied  where  ti  ;  detoua  >  on  fades.  bince  the  inner  pftion  is 
relatively  more  irratj'c,  t  is  b  ar  that  a  tv  «t  e»  dy  state  sus¬ 
tained  detonation  o*.  cun  s  tor  on  l.j  a  snot  dist  ice  and  his  deto¬ 
nation  is  in  tact  •.  rtnsieuf  phe  >■  n<u  on  . 
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EFFECT  OF  WEB  THICKNESS 


Two  samples  of  AAB-3189  propellant  were  tested  that  measured  8 -in. 

OD  by  1.5 -in.  ID  by  32 -in.  long.  These  wer:  fired  under  the  same 
conditions  and  with  the  same  type  of  instrumentation  as  that  shown 
in  Figure  IV-1.  The  data  from  the  two  tests  were  self -consistent. 
Figure  XV-8  shows  the  velocity -vs-distance  record  of  these  tests 
and  illustrates  the  behavior  of  detonation  velocity  in  a  hollow 
cylinder  that  has  a  web  thickness  (3.25  in.)  that  is  much  larger 
than  the  predicted  pscudocri  ti cal  web  thickness  (1.3  in .  ) . 

Throughout  the  entire  sample  length  the  velocity  of  the  wave  over 
the  outer  half  of  the  web  remains  constant  (4.8  tnm/usec)  .  At  the 
3/4-web  depth,  the.  same  velocity  is  maintained  for  10  in.  and  then 
the  velocity  increases  very  rapidly  to  reach  an  eventual  maximum 
of  9.5  mm/usec.  The  jet  velocity  increases  to  between  9.6  and 
9.5  mm/usec  after  12.5  in. 

Detonation  of  propellant  a  9.5  mm/psec  is  difficult  to  accept. 

Since  this  velocity  agrees  with  the  jet  velocity,  an  alternative 
explanation  for  the  high  velocity  registered  by  the  3/4-web  probes 
can  be  postulated.  Assume  that  the  jet  front  causes  ignition  of 
propellant  at  the  inner  surface  as  it;  proceeds  down  the  perfora¬ 
tion.  This  is  reasonable  because  of  the  high  temperature  and  high 
pressure  characteristics  of  the  jet.  Assume  further  that  because 
of  the  high  pressure  in  the  perforation  behind  the  jet  front  the 
propellant  will  burn  at  a  fast  rate.  It  follows  that  at  a  suffi¬ 
cient  distance  down  the  sample  enough  time  would  elapse  between 
the  passage  of  the  high-velocity  igniting  jet  front  and  the  arrival 
of  the  lower-velocity  reactive  wave  in  the  propellant  to  permit 
the  burning  propellant  surface  to  reach  the  inner  probes  before 
the  reactive  wave  in  the  propellant.  Since  the  probes  are  trig¬ 
gered  by  shorting  through  a  highly  ionized  medium,  they  could  be 
triggered  when  the  burning  surface  reaches  them.  The  apparent 
velocity  deduced  from  the  inner  (3/4-web)  probes  from  that  point 
on  should,  be  identical  to  the  jet  velocity,  because  the  vector 
describing  the  regression  of  the  burning  inner  surface  in  the  longi¬ 
tudinal  direction  is  equal  to  the  velocity  vector  of  the  jot  in 
the  same  direction  (assuming  a  constant  radial  burning -rate  vec- 
l  or  ) . 


5 See  paragraphs  entitled  Conclusions  £<u  u  discussion  of  tie 
critical  geometry  of  hollow  cylinders ,  and  the  definition  of 
p  s  eudoe  r i t  i c a 1 . 
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wlu'.ri'  Vi  i  ft  the  burn  in g  r;ite  (ips)  at;  pressure  P  (psig).  No  mea¬ 
surement  i.i  i  the  .static  pressure  within  a  cylindrical  jet:  eavi  ty 
exists,  so  there  is  no  way  Lo  determine  the  probability  that  such 
high  burning  rates  are.  possible. 


OONlVUiolONS 

Based  on  the  reported  findings  of  the  jetting  phenomena  test  series, 
the  following  conclusions  are  made.  By  all  considerations,  these 
conclusions  are  subject  to  revision  and  correction  pending  further 
study  of  this  complex  problem.. 

a.  Core  diameter  does  not  affect  the  abnormal  behavior  of 
the  detonation  wave  in  hollow  cylindrical  samples. 

b.  The  abnormal  behavior  of  the  detonation  wave  is  not 
caused  by  an  end-effect  mechanism. 

c.  The  principal  cause  for  abnormal  behavior'  of  the  deto¬ 
nation  process  in  hollow  cylindrical  samples  is  the 
web  thickness;  specifically,  the  size  of  the  web  that 
exceeds  a  pseudocrit ical.  value,  which  is  defined  a; 
that  size  below  which  no  transitory  sustainment  occurs. 

d.  The  mechanism  that  causes  abnormal  quenching  ol  the 
detonation  process  consists  of  a  radial  burning, 
directed  outward  from  the  inner  surf ace , which  pro- 
coeds  at  a  high  rate  and  reduces  the  web  size  to 
below  the  pseudocritieal  value. 

The  last  conclusion  expresses  the  present  interpretation  of  the 
test  results  which  does  nut  consider  any  other  perturbatio.  to 
steady-state  behavior  than  the  high  burning  rate  process.  If 
buch  a  process  can  continue  unabated,  clearly  no  classically 
supercritical  hollow  cylinder  is  possible.  However,  the  duration 
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Tim  '‘ill  ire  jo  I  ting  pro!)  I /mu  ii;  nut  well  oiiuigb  under:;  food  mi  i  !  s 
has  i  o  fundamentals  to  allow  further  .speculation  about  tills  and 
uiriur  pn.os  i  It  i  1  l  ti  cs.  Tlie.  primary  t|u««  1. i  mi  that:  remains  unanswered 
is  hov/  a  ;jet  traveling  at  about  twice  the  velocity  of  a  do  ton  at  ion 

noon  I  ion  in  the  web  of  a  hollow  sample  ram  I  i  titles  to  receive.  energy 

from  i  bat.  reaction  despi  l  e  the  continuously  I  tie  re  as  in#  scpara  t  ion 
of  the  two  front:;.  If  the  jet  velocity  exceeds  the.  detonation 
wave  velocity,  steady-state,  conditions  would  be  imposs  ibl  u  Ln 
hoi low -core  samples.  This  tael  atone  ru  Les  out  discussion  oi 

i  nil  ical  nf'/e  and  superci’i  I  leal  si'/.e,  when  speak  i  ng  of  the  liol  I  ow 

eyl  i  ncl'T.  Those,  terms  need  t.o  lm  mod  i  I.  j  od  by  us  inj<  tin1  prefix 
pseudo.  Thus,  pso.ndoori  t  i.r.nl  geometry  is  defined  ns  the  minimum 
geometry  in  which  detonation  can  be  sustained  for  a  minimal  dis¬ 
tance,  and  pseudosupei  critical  geometry  is  defined  a.s  any  geo¬ 
metry  larger  than  tiie  pseudocrit  ical  geometry. 


IV  I  . 


I  ge  L  i.d  -Vl  Op'-  !  !  an  i .  ll<»i?;i.ors  Hm.i  » 1  ns  j  ye.  i  |,-t  It? 

n  i  *  it  ly  Program  <  i*ri>;|4-c-:  ;-;c  si’ny  5 ,  '  i ' -  i  j  •  1 1  *.'/'* !  j  >**<••••<  >jitor»j 

Report  Ai'RIM.-TR-bS-:?,!  i  ,  /'ornjct  fp'niovt! 

O8f»t>~0  1(0*  )FI'f  24  November  5.44  j. 

IV--?.  Li.  win  ,  O.  R .  ,  P.  K.  '  zmnn  ,  and  W.  H.  Andersen, 

"Deflagration  Characteristics  of  Ammonium  Per  c  til  o 
rate  at  High  Pressures,"  Njnth  Symposium  ( Inter¬ 
national  )  on  Combustion,  Academic  Press,  New  Yorlc, 

piT358-3r,r'.  L%3.  ~ 


0.005 

0.004 


0.002 

0.003 


0 

.006 

0 

.006 

0 

.004 

0 

.003 

0 

.004 

Go 

4.26 

Go 

I  4.28 

i!u 

go 

Go 

4,2  5 

Go 

4.33 

Go 

4.26 

No 

go 

Go 

4.26 

Go 

4.31 

Go 

4.25 

Go 

4.30 

Go 

4.25 

Go 

'.29 

No 

go 

No 

go 

No 

go 

Go 

4.38 

Go 

4.33 

Go 

4.34 

Go 

4.32 

Go 

4.33 

Go 

4.37 

Go 

4.35 

Go 

4,36 

Go 

4.32 

Go 

4.37 

Go 

4.38 

Go 

4.41 

Go 

4.38 

Go 

4.36 

Go 

4.40 

Go 

4.44 

Go 

4.30 

Go 

4.41 

Go 

(No  Record) 

Go 

4.2  7 

Co 

4.39 

Go 

4.42 

Go 

4.41 

Go 

4,41 

Go 

4.38 

Go 

4,39 

Go 

4.38 

Go 

4.38 

!. ! .02 
! . I . 40 
!.i.35 
L.37 
M.38 
1.1.61 
l  .63 
1.1.64 
.1.42 
.  1 . 60 
.1 .31 
.1.36 
.1.33 
.1.32 
.1.34 
.1.24 
.1.47 
.1.52 


3.2.1.46 

3.2.1.59 

3.2.1.58 

3.2.1.69 
3.2.1.68 
3.2.1.45 
3.2.1.28 
3.2. 

3.2.1.71 

3.2.1 .72 

3.2.1.70 
3.2.1 .44 
3.2.1.27 
3  •>  .1.67 
3.2,1.54 
3.2.1.51 
3.2.1.50 
3.2.1.43 
3.2.1.66 

3.2.1.48 

3.2.1.49 


Average  Density  =  1.726  gm/cc 
Standard  Deviation  =  0.0048  gm/cc 


61 


r  J  S3  fC  K>  N  3 


:;s 

■;! 

1 


3 


f 


'  i ' ) }  >  1  ■,  ■.  V  .11.  I  j  I  <_;li  (,;y  I  iudi'i;:.. 


Mi  ,'m 

0  1  .  rr  i.i  l :  1 1  'i 

Or  V  ,  ;i  i  ,  r,n 
(.  ill.  ) 

l)<  *l  l : )  !  1  V 
(.  ) 

H  i  v ;  u  1  1 
cl 

1 

Avi'l  ago 

Uo  1  01  lo  1  i'lll 
Vi-1  or  i  t  y 
(  utiti/ 1 1  f ;  < 1  <  ■  ) 

S  t  ain  l;i  i  •<  i. 
l)<  •  \r )  i  •  i  i  •  >» i 
f)\ 

V.-  j  m  it  v 
( mm/ 1 1  s  <  *«) 

'iVs  t 

No  .  'J .  2  .  1 

.'ci 

0  .()():• 

1  .73  1 

!  N<>  go 

-  r  74 

■' .  49 

O  .(,)()  1 

1  .  /  7  H 

No  t’o 

74 

3 . 4  ) 

0  .002 

1  .  77  0 

No  go 

-  _ 

77 

2 .  49 

0 .007 

1.  .  72  2 

Mo  go 

-- 

7  3 

0  .007 

1  .  72.i> 

No  i.o 

7  4 

')  ^  1 J 1 , 

0 .008 

l. .  /;"■ 

•Jo  ;M1 

— 

84 

v  „  5  > 

0  .00  1 

1  .  77  4 

No  go 

HI 

0 .007, 

1  .  77  3 

No  \>ji 

84 

,  % ‘ » 

0  007 

1  .  77  7 

Mo  j,; o 

— 

H  ‘i 

.  v» 

0.001 

1  .  722 

No  j’.o 

.... 

-  - 

80 

.  55 

0  .002 

1,  .  72  7 

No  y\n 

-- 

78 

3  .  55 

0 .002 

I  .  721 

No  jo) 

— 

8  3 

2 .54 

0.001 

1  .771 

No  go 

-  - 

8  2 

■2  .  55 

0.007 

1  .  77. ! 

No  go 

— 

87 

7  .  I>  1 

0.002 

1.727 

No  go 

-- 

— 

94 

2-'.J 

0.0  01, 

1  7  O  C. 

No  o 

*  - 

94 

2  .  <»  1 

0.001 

1.720 

No  go 

--- 

— 

9.3 

2  .(.1 

0 .001 

I  .72  0 

’lo  go 

88 

2  .4  \ 

0  .002 

!.  .  724 

Mo  go 

97 

2  .<»l 

0.002 

1  .  774 

No  go 

-- 

94 

2  .cl 

0.003 

1  .724 

No  go 

91 

?  ,i,l 

0.002 

1  .  724 

Mo  go 

92 

2.nl 

O.OOj 

1  .72  3 

NO 

-- 

89 

2  .ill 

0.003 

1.721 

No  go 

. 

90 

2  -i>2 

0,003 

1  .  7  2  H 

N  o  go 

-- 

99 

2.1,2 

0.002 

1  .77  7 

No  go 

... 

100 

2.(,; 

0.002 

1,72/ 

No  go 

-- 

104 

2  .1,7 

o  .oo:' 

1  .  7  2  0 

No  go 

-- 

105 

:  .i)7 

0  .002 

1  .  72  0 

No  go 

-  -• 

1  01 

2  .i>7 

0.002 

L  .  722 

No  go 

LOO 

2  .  l,H 

0  .001 

1  .  778 

Go' 

4.12 

0.23 

1  03 

2  .  48 

o  .001 

1.727 

Oo 

4  .  L  9 

0.11 

102 

2  ,  OH 

0  .002 

1.7)8 

No  go 

-  - 

-  - 

98 

2  .  74 

,0 .004 

1  .7  33 

No  g;o 

-  * 

-  - 

113 

2 . 74 

0.091 

1  .729 

.  <  if) 

4.21 

0.  It) 

1  11 

2.74 

0  .01)1 

1  .739 

Oo 

4.25 

9 . 1  0 

114 

2 . 74 

0 .00  1 

1  .  724 

(  in 

4.25 

0.10 

1  1.4 

2.74 

0.00  3 

1  .728 

No  go 

1  1.2 

2 .74 

0.00  L 

1  .  /:■/ 

Gu 

4.2  5 

0 .03 

1  09 

2 . 74 

0.001 

1  .72/ 

No  go 

1 13 

2 .74 

0.001 

1.727 

(in 

4.;  9 

0 . 04 

1  08 

2 .74 

0 . 00  L 

l  .  7  24 

i » o 

4.2  4 

0.99 

1  1  0 

2 . 74 

0 . 00  1 

1  .  724 

IjO 

l  no  record ) 

-  -  1 

107 

2 . 79 

0  .002 

1.731 

G.) 

4 . 2  J 

0.11 

120 

2.79 

0 .004 

1.720 

Go 

4  .29 

0. 1  3 

i  I  9 

2  .80 

0.002 

1  .  7  30 

tin 

h.  7  7 

0.04 

1  l.H 

2  .80 

0.003 

1  .729 

Go 

4.2. 

0  .  L9 

1  21 

0 .  Oil  8 

1  .72  3 

Go 

4 . 1.  4 

0.10 

“117 

Avcnm'.c  Density  =  1./24  gin  4  r 

!  r-vi.it  i  on  -  O'.OOJI  ^ 


i 

\ 


Ni  Si  kj  r  j 


Avuras;*  Density  =  J. .  725  gro/cc 
Standard  Deviation  =  0.0023  gra/<  c 


63 


9 

J  P.  f  .. 

.'{Ms’ 

Si  am  i . 1 1  <  I 

Mralt  l)('V  i;i  I  J  «’ ? 1 1  Ui'iis  J  t  y 


in,  ) 

(.  in.) 

(  yt\u/<  ■< 

4 . '/ ') 

O.OOI 

1  .  /  37 

h  .  ,1 S 

0  .02  7 

i  .  7  •;  3 

4 .  a  r> 

0.04/ 

1  .7  30 

4.33 

0.040 

]  .  720 

4.3/ 

0.034 

1  .  /  3  1 

4 . 3K 

0.030 

1.720 

4.4/ 

0  043 

1  .  7  30 

4 . 48 

0.0  30 

1  .  7  32 

4 . 44 

0.038 

1  .  7  30 

4  .  r,o 

0.040 

1  .  7  34 

4  .  ‘>0 

0.040 

1  .  73  1 

4 . 40 

0  .0  44 

1.731 

4.  44 

0 .038 

!  .  720 

4 44 

0.032 

1  .  728 

4.40 

0.044 

1.731 

4.44 

0.04] 

1  .  730 

4 .  (i4 

0.02.4 

1.  728 

4.07 

0.054 

1  .  7  30 

4.4/ 

0 . 043 

L  .  727 

4.  /!) 

0.0  70 

1  .  720 

4.74 

0  052 

1  .  ’ 

i  {  i-  'i  i 1 ,  r  "j  i  <.  i  1 , ;  \  5 . *  1  ■ )  i  1 J 1  *  «  ,it  i  *1  ^  - 


Av(»  »  ■ 

: :  1 1  }.  J  : ;  l)>  -  I  mu  !  i  <  >n 


o  A. 

T-oH 

V<  >  ]  1 1 1  i  1  y 

(  MIII(/||JilM  '  ) 

'IV.'i  i  Nil, 

No  jui 

. 

3  ..  3.2.  2 

No  go 

3 . 3 . 7 . 3 

No 

-  • 

3  .  3 . 7  ..  8 

Nil  (i,  o 

3.3.  ,2 . 0 

No  go 

1.3.2. 10 

No  go 

3  .  i .  7  .  1  I 

No  go 

- 

3 .3.2.17 

No  t'o 

3.3.2.  i  3 

No  j*  i » 

. 

i  .  3 . 7 . 4 

No 

4.34 

3 . 3 . 7  .  J, 

4.7.4 

3 .3.2  ..  4 

No  go 

3.3.2 .24 

Co 

4.32 

3.3  2.4 

Go 

4.74 

3.3.2.10 

Go 

4.  30 

3.3.2./ 

Go 

4.20 

3.3.2.14 

Go 

4.34 

3.3.2.21, 

Go 

4.31 

3. 3.2 .20 

flu 

4.32 

3 , 3 . 2  J  8 

0  o 

No  data 

3.3.2..1  0 

Go 

4.4,1 

3.3.2.17 

Table  V-V.  Batch  ^Ell-110  1. /5-in. -Thick  Rectangles* 


(icwv  T  ry 
f  III.  ) 

Width 
f  in.) 

Result s 

Avr  rage 

Dp  t  oria  t  i  on 

V  i- 1  o<  i  l  y 
(  mm/' g  sot  ) 

Test  N  o . 

2 . 40 

5 .  CO 

No  go 

_ 

3.3.2.78 

2  40 

5.75* 

No  go 

3.3.2.77 

2  . 49 

5.75 

No  go 

3 . 3 . 2 „ 49 

2  .40 

5.875 

No  go 

3.3.2.  50 

2  .  46 

5.8  75 

No  go 

-- 

3  3.2.51. 

2.7* 

4.00 

No  "o 

3.3.2.48 

2  7! 

4.00 

No  go 

3.3.2.42 

2  7  , 

4.00 

Go 

4 .  i  9 

3.3.2.37 

2  7  3 

4.125 

No  go 

3.3.2- 44 

2  7  3 

4.1  ?,5 

No  go 

-- 

3.3.2.47 

2  74 

0.25 

No  go 

3.3.2.45 

2.70 

4.50 

Go 

4 . 38 

3 .3.2.44 

2  70 

7  .00 

Go 

4 . 24 

3 . 3 . 2 . 43 

2  So 

7.75 

Go 

4.29 

3. 3.2 .41 

2  03 

9  .00 

Go 

4 . 29 

3.3.2.38 

2  96 

0 . 50 

Go 

4.13 

3.3.2. 40 

2  09 

10.00 

Go 

4 . 24 

5.3. 2.34 

3  .  U5 

11.75 

Go 

4,31 

3.3.2.35 

a 

i 


s 

? 

9 

l 

I 


! 

» 

f 

i 


Sample  length  -  /3-5  in 

Aye  ra,e,e  densf  t  y  -  1.72a  yin  y< m r 

6h 


<)U 


4.  2 5 
4 .  2  5 
4.  I  2 
4 .  i  2 
4 . 0() 
4 . 00 
3.88 
3.88 
8 ,  73 

3.75 
3.02 
3.82 

5..  75 

5.75 
5.62 
5.62 
5.50 
5.50 
5.38 
5.38 
5.25 
5.25 
5.25 

r- 

**  m  Hm 

5.12 


l  >:  i  ji1;  i  I 

U  1  H  I  Kill 
t  i  1 1. ,  ) 


.13 
17 
2 . 5 
4 
0 
0 
0 
( I 

0 

o 

0 

0 


8.5 

4.5 
5 

0 

0 

0 

0 

0 

0 

0 

K* 

0 


r :  ic  f« <;;co!'  is ;  Mjcriioiui 


i'ni'  i  n  i  N! ; 


vi!i)!  curri'ENT 


.i„  L.  Mi  i  ;m  I  ,  tliiomma.t  .Special  is » 
Aero  iet  -( Leri  e  ra  )  Corporal  f. on 
Sacrament  o ,  Cal  l"  t  orn  i a 


INTPoniTCT  ION 

This  report  is  on  that  phase  of  the  SOPHY  program  concerned  with 
the  evaluation  or  microscopic  methods  for  determining  binder  void 
content.  The  objectives  of  this  study  were:  (1)  to  determine 
whether  small  voids  in  the  binder  could  be  recognized  with  the. 
microscope. .  (?)  to  evaluate  microscopic  methods  for  determining 
loL.nl  hinder  void  volume,  in  a  propellant  sample,  and  (3)  to  deter¬ 
mine  the  size  range  oi  voids  in  a  propellant  sample  with  the  micro¬ 
scope  . 


VOID  RECOGNITION 

In  Objective  1,  the  small. st  void  previously  observed  was  100  g. 
Production  samples  contain  voids  of  a  minimum  of  15 p  in  diameter. 

A  difference  in  the  ability  to  recognize  voids  in  this  smaller 
size  range  results  from  the  method  of  sample  preparation.  The 
propellant  is  microtomed  into  thin  sections  of  ;  predetermined 
thickness.  If  this  thickness  is  30  g,  any  void  larger  than  30  u 
in  diameter  will  pasc  completely  through  the  thin  section  and  is 
easily  seen  in  transmitted  light  as  clear  white  circles,  as  shown 
in  Figure  VI -1.  This  section  is  dry  mounted  for  demonstration 
purposes;  normally  the  section  must  be  oil  mounted,  as  in  Fig- 
u  s  VI-2  b,  c,  and  d.  If  a  void  smaller  than  30  p  is  present,  it 
will  be  totally  included  within  the  section.  If  oil  does  not 
penetrate  the  void  it,  in  theory,  will  appear  as  a  black  opaque 
spot  such  as  the  dry  void  inside  the  AP  crystal  in  Figure  VI-2  c. 
The  void,  is  nearly  opaque  because  it  acts  as  a  spherical  mirror 
in  transmitted  parallel  light.  If  parallel  light  is  made  to  con¬ 
verge  on  the  void  with  the  swing -in  substage  condenser,  a  white 
spot  appears  in  the  center  of  the  void  as  shown  in  Figure  VI -2  d. 

In  this  way  the  void  can  be  distinguishe^a  from  spherical,  com¬ 
pletely  o pa sue  aluminum  particles.  This  optical  interpretation 
is  confirmed  when  the  same  void  is  examined  after  the  oil  has  pene¬ 
trated  the  void  and  it  is  seen  to  transmit  light.  This  reasoning 
will  apply  on  voids  down  to  about  1 p,  in  diameter. 
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There  i.r;  .*  unique  morpho  1  <>g  i  <:a  L  nt.na-  turu  in  the.  |iro|ii>.lJ.ant  around 
these  ;  mal L  voids  shown  in  Figure  VI -2  b  and  n.  The  void  is  sur- 
t^ivi i nrl  a  *■* i  i  r\f  ma  ~ap  ‘/hich  hs  iLti  turn  wiirr'O'iindLcd.  by  ci  ch.c. LI, 

of  cl  ar  binder.  Two  ol  these  structures  are  shown  in  Figure 
VI -2  b.  In  Figure  VI-2  c,  the.  void  is  missing  and  a  bail  of  MA- 
i‘  is  off  center  in  the  shell  of  binder.  This  structure  is  usc- 
tal  for  finding  the  small  voidr  when  counting. 


THEORETICAL  CONSIDERATIONS  OF  VOID  VOLUME  AND  SIZE 

DISTRIBUTION 

The  determination  of  u  volume  from  microscopic  study  of  a  two- 
dimensional  field  was  mathematically  developed  by  Roaiwal  in  1898. 
Numerous  tests  were  conducted  and  a  review  of  the  procedure  by 
E.  S.  Larsen  and  F.  S.  Miller  was  published  in  the  Journal  of 
Mineralogical  Society  of  America,  Vol.  20,  No.  4,  April  1935,  lti 
which  they  state  "that  the  linear  method  of  Rosiwal  in  easily  cap¬ 
able  of  an  accuracy  within  1  or  2  percent  for  each  constituent." 

In  SOPHY  propellant,  the  voids  are  considered  to  consist  of  a 
series  of  different  size  spheres,  randomly  distributed  with  a 
loose  packing.  This  condition  is  slightly  different  from  the 
Rosiwal  problem  tor  an  intimately  packed  crystalline  mosaic; 
however,  the  technique  should  be  valid. 

The  distribution  of  void  sizes  was  given  a  partial  mathematical 
development  by  R.  Farris  of  Aerojet -General  by  limiting  the  prob¬ 
lem  to  spherical  voids  of  one  size,  and  is  as  follows: 
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(VI-1) 


If  a  line  of  normal  incidence  intersects  the  sphere  at  distance  of 
radius  r,  the  circle  formed  by  this  radius  represents  the  locus 
of  all  chords  of  length  P.  So  that 


when  P  =  2  R,  r  =  0 
and  P  =  0 ,  r  =  R 


(VI -2) 


If  the  probability  of  intersecting  any  point  on  the  circle  is  the 
same,  then  the  probability  P  of  measuring  a  chord  length  from 
t,  -  0  to  P  is 


P<  P)  -P(0) 


(VI-3) 


where  D  is  the  diameter  of  the  void.  Therefore,  Farris  concludes, 
if  a  histogram  is  plotted  on  log -log  paper,  it  should  have  a  slope 
of  2  and  an  intercept  of  1/d2. 

A  graphical  model  was  constructed  with  spheres  of  one  size.  The 
Rosiwal  equation  100  x  /(x  +  y  +  z)  was  used  to  calculate  the 
volume  of  voids  and  a  plot  was  made  to  determine  the  slope.  The 
results  were  slightly  in  error,  presumably  because  of  the  small 
site  of  the  graphical  model,  but  sufficiently  accurate  to  proceed 
with  this  experiment. 


VOID  VOLUME 

An  experimental  run  was  made  on  Sample  82,  which  is  visually  simi¬ 
lar  to  the  sample  in  Figure  VI -2  a,  to  determine  the  void  volume. 
Two  different  microscopic  methods  were  used.  Both  made  use  of 
the  linear  method  of  traversing  the  thin  section  and  recording  the 
lineal  intercepts  of  propellant  and  void  along  a  line.  In  one 
method,  an  eyepiece  micrometer  was  used  to  measure  each  intercept 
distance  and  resulted  in  a  tabulation  of  each  void.  Iu  the  other 
method,  the  Leitz  Integrating  Stage  was  used  which  gives  a  cumula¬ 
tive  reading  of  all  intercepts  ana  the  number  of  voids  was  not 
recorded. 
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On  the  test  sample,  the  length  of  the  linear  traverse  was  7.4“  cm 
and  flic  percent  void  volume,  calculated  was  i  1  ,k%  using  i  iie  eyep  ieue 
iriicromecer.  A  short  traverse  of  1.3  cm  gave  11.0%  void  volume 
Using  the.  integrating  stage.  A  density  detenu  ina  L  i  on  on  the  hulk 
snMiples  by  an  oil.  displacement  method  on  an  analytical  balance  gave 
a  void  volume  of  12. 0%.  The  1%  discrepancy  should  be  due  to  a  lew 
larger  pores  in  the  bulk  sample  that  were  deliberately  avoided  in 
microtoming  the  thin  section.  The  method  using  the  integrating 
stage  required  20  min  while  1  he  eyepiece  micrometer  method  re¬ 
quired  all  day. 


VOID  SIZE  DISTRIBUTION 

When  the  above  eyepiece  measurements  of  individual  'mids  were 
plotted  on  log-log  paper,  the  slope  of  the  line  was  approximately 
1,  which  indicates  a  distribution  of  void  sizes.  The  mathematical 
procedure  developed  by  R.  Farris  could  be  extended  to  determine 
the  void  size  distribution  from  this  single  line,  but  would  require 
several  weeks.  Lacking  this  method,  an  experimental  method  was 
attempted.  In  studying  the.  voids  at  high  magnification,  it  was 
noticed  that  the  walls  of  the  voids  passing  through  the  thin  sec¬ 
tion  were  cusped  up  or  down,  and  a  few  were  vertical.  From  the 
geometry  of  these  voids,  it  is  apparent  that  cusped  walls  are 
formed  by  minor  circles  in  a  sphere  and  the  vertical  walls  are  from 
a  great  circle.  Thus,  by  scanning  the  thin  section  and  measuring 
the  diameter  of  holes  with  verticaL  walls,  the  void  size  distribu¬ 
tion  can  be  recorded.  Voids  that  do  lot  pass  through  the  section 
always  give  a  great  circle  projection  and  all  of  these  can  be 
measured.  To  find  sufficient  voids  for  measurement,  several  thin 
sections  were  cut  from  different  parts  of  the  block  sample  No.  82. 

A  total  of  207  voids  were  recorded  with  a  size  range  from  1204 
to  17  (A.  The  data  were  statistically  reduced  and  placed  on  normal 
probability  paper,  Figure  VI-3.  There  is  a  dog  leg  in  the  curve 
at  the  median  point  and  another  where  three  very  large  voids  appear 
Two  dog  legs  were  expected  on  the  basis  that  there  is  a  bias  in 
counting  a  greater  number  of  voids  that  do  not  pass  through  the 
thin  section  and,  secondly,  as  pi'eviously  reported,  some  of  the 
very  large  voids  were  obviously  the  result  of  coalescence  of  the 
two  smaller  voids.  Interestingly,  the  bias  on  the  graph  of  the 
smaller  voids  is  in  the  wrong  direction  and  the  size  is  about  twice 
the  thickness  of  the  thin  section.  This  corresponds  to  the  size 
range  of  the  voids  which  are  enclosed  in  the  unique  binder  shells 
shown  in  Figures  VI -2  b  and  c.  Disregarding  the  few  cases  of 
coalescence,  it  appears  that  there  mav  be  a  bimodal  distribution 
of  void  sizes. 
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Figure  VI-3.  Distribution  of  Pore  Sizes. 
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. .  ShockWave  Attenuation  in  12-in.  Squares  of 
Stacked  PMMA  Plates ,  TNT  Booster . 
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Figure  VII-5,  Shock  Wa-  e  Attenuation  in  24-in.  Squares  of 
Stacked  PMMA  Plates,  TNT  Booster. 
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Figure  VII  -6 


Shock  Wave  Attenuation  in  42-in. 
Slacked  PMMA  Plates,  TNT  Booster. 
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Table  VIII -I „  Circular  Cylindrical  Acceptor. 


Test 

No. 


3. 3. 7. 3 


3. 3. 7. 5 


Time 
( |j,sec) 


5.0 
10.0 
15.0 
20.0 
25.  0 
30.0 
35.0 


5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 


5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 


5.0 

10. 

15. 

20. 

25. 

30.0 

35.0 


5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 


h 

(  i  n .  ) 


.34 
.30 
.51 
.07 
.42 
.48 
1.80 


0.06 

-0.71 

•0.20 

0.81 

1.03 

1.12 

0.83 


-0.83 
-0.57 
-0.  11 
0.33 
0.57 
0.74 
1.45 


(  i*v . ) 


2.93 

3.34 

4.01 

4.48 


2.75 

3.15 

3.65 

4.47 


3.12 

3.60 

4.03 


2.08 

2.79 

3.25 
3.68 

4.26 
4.91 
5.  10 


X1 

( in . ) 


0 

0.67 
1.24 
1.85 
2.50 
3.23 
3.  95 


4.14 

5.13 


3.54 

4.14 


0.10 

0.93 

1.30 

2.07 

2.67 

3.20 

3.94 


0.00 

0.52 

1.26 

2.46 

3.59 

4.63 

5.57 


Average 

Shock 

Pressure 

Wave  Area 

( kbar ) 

( in. )2 

31.43 

11.1 

16.04 

15.85 

9.91 

22.60 

6.35 

25.60 

3.97 

24.86 

2.26 

18.99 

1.03 

10.94 

52.79 

11.40 

33.27 

15.49 

25.70 

24.90 

21.90 

28.50 

20  ,25 

3 ‘.15 

28  .00 

33.99 

34.30 

32.45 

41.33 

4.29 

29.74 

26.38 

20.77 

29.80 

14.59 

21.76 

10.50 

24.90 

7.40 

28.30 

5.55 

32  ,47 

34.54 

9.24 

24.70 

10.82 

19.40 

19.50 

15.66 

22.68 

13.01 

26.76 

11.09 

33.31 

10.99 

33.91 

55.10 

16.34 

36.05 

29.87 

29.53 

38,43 

34.35 

35.80 

41.86 

33.07 

50.58 

31.56 

60.80 

31.27 

'uur*Li 


Tabic  V1II-I.  Circular  Cylindrical  Acceptor  (G<«nt. „  ) 


Test 

No. 

Time 
( psec ) 

h 

(in.) 

(in.) 

V 

“1 

(  in  .  ) 

Average 
Pressure 
( kbar ) 

Shock 
Wave  Area 

( in  .  ) ^ 

3. 3. 7. 6 

10.0 

0.30 

1.94 

1.13 

" 

28.02 

13.49 

15.0 

0.50 

2.52 

1,72 

20.96 

20,81 

20.0 

0.97 

2.80 

2.41 

17.23 

24.01 

25.0 

1.37 

3.11 

3.03 

15.34 

28.37 

30.0 

1.82 

3.34 

3.63 

14.59 

32.09 

35.0 

2.33 

3.50 

4.22 

14.70 

35.35 

1 

5.0 

0.51 

0.75 

0.77 

47.68 

2.31 

j  1 

10.0 

0.44 

1.71 

1.09 

38.20 

11.39 

II 

15.0 

0.68 

2.35 

1.64 

32.17 

20.57 

II 

20.0 

1.56 

2.39 

2.59 

28.06 

20.48 

HSfirauH 

25.0 

2.42 

2.42 

3.50 

25.53 

20.42 

II 

30.0 

3.20 

2.47 

4.32 

24.18 

20.93 

■1 

35.0 

3.84 

2.56 

5.01 

23.62 

22.32 

3. 3. 7. 8 

10.0 

-0.17 

1.98 

0.76 

21.99 

13.05 

15.0 

-0.02 

2.56 

1.31 

16.60 

19.80 

20.0 

0.48 

2.76 

2.03 

13.03 

20.96 

25.0 

0.94 

2.92 

2.69 

10.62 

21.50 

30.0 

1.47 

2.99 

3.36 

8.79 

20.68 

35.0 

2.04 

2.99 

4.02 

7.40 

19.00 

3.3,7.12 

10.0 

0.  56 

1.79 

1.45 

'  .  v 

10.11 

15.0 

0.50 

2.53 

1.87 

' 

18.40 

20.0 

0.74 

3.00 

2.39 

■QHH 

25.48 

25.0 

0.89 

3.48 

2.70 

16.39 

36.42 

3.3.7.14 

5.0 

0.46 

0.84 

0.73 

58.11 

1.36 

0.93 

1.81 

20.88 

■ms 

■atB 

0.23 

2.78 

1.59 

19.77 

0.69 

3.10 

3.35 

15.27 

28.14 

1.42 

3.  16 

3.19 

13.11 

27.54 

30.0 

1.72 

3,50 

3.71 

12.57 

33.20 

35.0 

2.31 

3.59 

4.33 

13.02 

35.44 

Table  VII1-II.  Rectangular  Accentor 


Teat 

No. 


3-3. 8. 2 


Time 
( nsec) 


5.0 
10. 0 
15.0 
20.0 
25.0 
30.0 
35.0 
40.0 
45.0 
50.0 
55.0 
60.0 
65.0 
70.0 


3. 3.8.3  5.0 


h 

(in.) 


R 

(in.) 


30.0 

35.0 

40.0 

45.0 

50.0 

55.0 

60.0 

65.0 

70,0 

75.0 


0.67 

1.44 

2.09 

2.81 

3.54 

4.03 

4.60 

5.24 

5.48 


-1.93 

-1.14 

-1.17 

0.01 

-3.87 

1.45 

0.38 

0.59 

-1.91 

4.19 

3.19 
2.55 


4. 

.16 

4. 

.16 

4. 

.  35 

4. 

,45 

4. 

,77 

4. 

,92 

5. 

05 

5. 

.31 

5. 

.48 

5. 

.66 

6, 

.06 

6. 

.38 

6. 

,66 

7. 

,29 

6, 

.32 

5, 

.47 

4. 

,03 

5. 

.85 

5. 

,47 

6. 

,85 

6. 

,52 

11, 

.27 

6. 

.  8i 

8, 

.75 

9, 

.42 

12. 

.08 

7. 

.61 

9, 

.42 

10, 

.96 

k 

( in. ) 


-1.69 

-1.80 

0.02 

0.61 

0.99 

1.76 

2.62 

3.52 

4.13 

3.88 

4.95 

5.48 

6.38 

7.55 


(  in.  ) 


2.91 
3.89 
3.01 
3.29 
3.78 
3.  90 
3.93 
3.93 
4.22 
5.33 
5.16 

5.52 

5.52 

5.52 


Average 
Pressure 
( kbar) 


Shock 
Wave  Area 
(in. )2 


4.89 
4.57 
3.94 
4.70 
3.48 
5.31 
4.20 
4.36 
4.51 
4.44  I  4.75 
5.30  4.80 

.08  4.88 
.73  5.14 

.18  5.52 

.81  5.77 


120.02 

22.9  2 

111.89 

32.84 

107.85 

40.49 

101.13 

41.81 

93.42 

40.00 

89.11 

39.23 

86.52 

38.73 

84.86 

38.02 

84.24 

37.28 

84.01 

36.36 

85.89 

35.88 

88.83 

35.28 

93.43 

33.84 

100.09 

34.76 

124.30 

31.00 

112.69 

36.88 

106.64 

37.48 

91.40 

36.40 

86.77 

37.36 

78.51 

35.04 

77.15 

35.92 

72.41 

34.00 

74.2  7 

35.44 

73.62 

34.32 

75.43 

34.08 

77.91 

30.80 

85.36 

34.60 

90.39 

33.80 

99.23 

33.44 

Figure  IX-2.  Flyer  Plate  Velocity  vs  Wt/Sq  In.  of  Det^ereet- 
(Aluminum  Flyer  Plate  1/8  in.  by  5  in.  by  5  in,; 
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Table  X-II.  PeaK  Side-On  Overpressure  and  Inorlse 


Table  X-VI.  Peak  Side-On  Ovemressure  and  Immalse 


'Computed  froat  shock  velocity  ncasr'rerp.eTits  using  Rankine-Hugoniot  equation. 
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Side-On  Overpressure  and  Impulse  (Concluded). 


Table  X-XII,  Peak  Side-On  Overpressure  and  Impulse  (Continued). 


Table  X-XIII.  Peak  Side-On  Overpressure  and  Impulse. 


Table  X-XIII.  Peak  Side-On  Overpressure  and  Imoulse  (Concluded) 


Table  X-XIV.  Peak  Side-On  Overpressure  and  Impulse. 
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Table  X-XIV.  Peak  Side-On  Overpressure  and  Impulse  (Concluded) 


Table  X-XV.  Peak  Side -On  Overpressure  and  Impulse 
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Table  X-XV.  Peak  Side-On  Overpressure  and  Impulse  (Continued). 


Table  X-XV.  Peak  Side-On  Overpressure  and  Impulse  (Concluded). 


.a 


..o  x> 

C  H  rl 

co  Q  O  O 

on  oC  o  ."t  rs.  tu  rj 
«  • 

0^4  *  o 
yO-jo'4<7'Z 


w  * 

o  je  H 

ss.s.f 

o  *>  -C  R 

•4H  p 

U  iJ  4J  4J  r-« 

.5  5  5*22 

lo  ‘H  N 
VMM 
W  «  <-» 

lliflii  V  W 
•»  O  0  o  0  «l  ** 
1  U  U  U  O  «  V 


« 

fed 


x 

'3 

> 


M 

41 

« 

5 


V 


&  3 


H 


114 


•  Vi  ■**««." 


Table  X-XVI.  Peak  Side-On  Overpressure  and  Impulse. 
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Table  X-XVII.  Peak  Side-On  Overpressure  and  Impulse  (Continued) 


XVIII.  Peak  Side-On  Overpressure  and  Impulse  (Concluded) . 


Table  X-XIX.  Peak  Side-On  Overpressure  and  Impulse. 
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Table  X-XIX.  Peak  Side-On  Overpressure  and  Impulse  (Continued). 


Table  X-XIX,  Peak  Side-On  Overpressure  and  Impulse  (Concluded) 


Peak  Side -On  Overpressure  and  Impuls 


Table  X-XX.  Peak  Side-On  Overpressure  and  I^pulsi.  ^Conclude 


Table  X.-XXI,  Peak  Side-On  Overpressure  and  ircul 


Table  X-XXI,  Pe  i  Side-On  Overpressure  and  Impulse  (Continue. 


mm 


Table  X-XXII.  Peak  Side-On  Overpressure  and  Impulse 


Teile  X-XXII.  Peak  oide-On  Overpressure  and  Itroulse  (Cent 
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Table  X-XXII.  Peak  Side-On  Overpressure  and  Impulse  (Concluded). 
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